Thermodynamics of the brane 
in Chern-Simons matter theories with flavor 

Niko Jokela.rl Javier Mas, 1 !] Alfonso V. Ramallo,^ and Dimitrios Zoakos^ 



1 Departamento de Fisica de Particulas 
Universidade de Santiago de Compostela 
and 

Instituto Galego de Fisica de Altas Enerxias (IGFAE) 
E- 157 82 Santiago de Compostela, Spain 

2 Centro de Fisica do Porto 
and 

Departamento de Fisica e Astronomia 
Faculdade de Ciencias da Universidade do Porto 
Rua do Campo Alegre 687, ^169-007 Porto, Portugal 



Abstract 

We study the holographic dual of flavors in a Chern-Simons matter theory at non- 
zero temperature, realized as D6-branes in the type IIA black hole dual in the ABJM 
background geometry. We consider both massive and massless flavors. The former 
are treated in the quenched approximation, whereas the massless ones are considered 
as dynamical objects and their backreaction on the geometry is included in the black 
hole background. We compute the holographically renormalized action of the probe 
by imposing several physical conditions. In the limit of massless flavors the free energy 
and entropy of the probe match non-trivially the first variation of these quantities 
for the backreacted background when the number of flavors is increased by one unit. 
We compute several thermodynamical functions for the system and analyze the meson 
melting phase transition between Minkowski and black hole embeddings. 
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1 Introduction 



Recent studies of Chern-Simons matter theories in three dimensions by holographic tech- 
niques have provided non-trivial examples of the AdS/CFT correspondence [UG] which 
could be of great help to shed light on the dynamics of some strongly coupled systems 
in condensed matter physics. The paradigmatic example of these systems is the Aharony- 
Bergman-Jafferis-Maldacena (ABJM) theory constructed in [3], based on the analysis of [HE], 
where the supersymmetric Chern-Simons matter theories were proposed as the low energy 
theories of multiple M2-branes. 

The ABJM theory is an M = 6 super Chern-Simons gauge theory in 2+1 dimensions with 
gauge group U(N)k x U(N)_k with opposite level numbers k and —k. In addition to the 
two gauge fields, this theory contains two pairs of chiral superfields which transform in the 
(N, N) and (N, N) bifundamental representation. When N and k are large the theory admits 
a geometric description in terms of an AdS^ x CP 3 with fluxes in type IIA supergravity which 
preserves 24 supersymmetries. The study of this theory and its generalizations has uncovered 
a very rich structure and has provided new precision tests of the AdS / CFT correspondence 
(see [6H9] for reviews of different aspects of the Chern-Simons matter theories). 

The ABJM theory can be generalized in several directions. In this paper we will consider 
the addition of fields transforming in the fundamental representations (N, 1) and (1, N) of the 
U(N) x U(N) gauge group. It was proposed in [TOlfTT] that these flavors can be incorporated 
in the holographic dual by considering D6-branes that fill the AdS^ space and wrap an RP 3 
submanifold of the internal CP 3 space. These configurations are M = 3 supersymmetric. 
When the number Nf of flavors is small one can adopt the so-called quenched approximation, 
in which the flavor D6-branes are considered as probes in the AdS± x CP 3 geometry. This 
approach has been followed in [T2T - H5] . 

In [12] a holographic dual of ABJM with unquenched flavor was found by considering a 
large number Nf of flavor D6-branes which are continuously distributed in the internal space 
in such a way that M = 1 supersymmetry is preserved. To find the unquenched solution 
one has to solve the equations of motion of supergravity with brane sources, which modify 
the Bianchi identities of the forms and the Einstein equations. If the branes are localized, 
the sources introduce Dirac 5-functions in the equations, which makes the problem very 
difficult to solve. For this reason we will follow the approach initiated in [T7] and study the 
backreaction induced by a smeared continuous distribution of flavor branes. This procedure 
has been successfully applied to add unquenched flavor in other holographic setups [T8l - f20] 
(see [2T] for a review and more references). As the smeared flavor branes are not coincident, 
the flavor symmetry for Nf branes is U(l) Nf rather than U(Nf). Moreover, since we are 
superimposing branes with different orientations in the internal space, the corresponding 
supergravity solutions are generically less supersymmetric than the ones with localized flavor. 
The unquenched solutions with smeared flavors are much simpler than the localized ones and, 
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in many cases the solutions are analytic. 

The unquenched solution of type IIA supergravity found in [TB] includes the backreaction 
effects due to massless flavors. The corresponding ten-dimensional geometry is of the form 
AdS± x A^6 ; where A^6 is a compact six-dimensional space whose metric is a squashed version 
of the unflavored Fubini-Study metric of CP 3 . In this solution the deformation introduced 
by the flavors is encoded in the squashing factors, which are constant and depend non- 
linearly on the number of flavors (although the sources of supergravity are linear in Nf). 
Notice that the backreacted metric contains an Anti-de Sitter factor. This is related to the 
fact that the dual Chern-Simons matter theory has conformal fixed points even when the 
flavors are added (see [22J for a verification of this property in perturbation theory). On 
the gravity side this conformal behavior is responsible for the regularity of the metric at the 
IR, contrary to other solutions with unquenched massless flavors [21]. It was checked in [16] 
that this solution captures rather well many of the effects due to loops of the fundamentals 
in several observables. In particular, it matches remarkably well with the behavior of the 
effective number of degrees of freedom of the flavored theory in the Veneziano limit, which 
was computed in the field theory side using localization in [23] . 

In sharp contrast to what happens to other flavored backgrounds obtained with the smear- 
ing method (see, for example, those of refs. pT8H20] ) . our supergravity solution has a good UV 
behavior and, since the metric has an Anti-de Sitter factor, we are dealing with a geometry 
for which the holographic methods are firmly established and it is possible to apply a whole 
battery of techniques to perform a clean analysis of the different flavor screening effects. In 
particular, as it is shown below, it is straightforward to add a further temperature defor- 
mation to the flavor deformation and to construct a black hole which contains the effects 
of massless flavors. This is simply done by including the standard blackening factor in the 
Anti-de Sitter part of the metric, without modifying the internal space M.§. We can then 
compute different thermodynamic quantities for this flavored black hole. 

When flavor branes are embedded in a black hole geometry the system undergoes a first 
order phase transition when the branes fall into the horizon [2H[25]. On the field theory 
side this phase transition corresponds to the melting of mesons in a deconfined plasma. The 
analysis of the influence of unquenched flavor in this melting transition is clearly a very 
interesting problem. However, in order to have a complete understanding of this problem in 
the holographic setup one has to find a black hole solution containing the full backreaction 
of massive flavors, which is very hard to find. In this paper we will adopt a more modest 
approach and consider a small number of massive flavors and a large number of massless 
quarks. The latter will be included in the background, while the massive fundamentals will 
be treated in the quenched approximation. Accordingly, we will consider a D6-brane probe 
in the non-zero temperature version of the background found in [16] and we will study its 
thermodynamic properties, following the same methodology as the one employed in [23] for 
the D3-D7 and D4-D6 systems. 
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The action that governs the dynamics of our D6-brane probes contains a contribution 
from the Dirac-Born-Infeld (DBI) and Wess-Zumino (WZ) terms. This probe action must 
be renormalized holographically in order to get finite answers for the different thermodynamic 
functions. At zero temperature one can adopt a gauge for the RR seven-form potential CV in 
which the two terms of the action cancel with each other on-shell for the kappa symmetric 
embeddings of the probe. At non-zero temperature the on-shell action of the probe in this 
gauge is finite, and the only freedom left by the holographic renormalization is the addition 
of finite counterterms. These finite terms can be fixed by imposing regularity of CV at the 
horizon and by requiring that all the thermodynamic functions for the probe vanish for 
infinitely massive flavors, as they can be integrated out. 

Once the action of the probe is fixed in this way, we should verify that it satisfies a non- 
trivial compatibility condition with the background. Indeed, let us consider a probe for a 
massless flavor. In this massless limit the quarks introduced by the probe are of the same 
type as those of the background. Thus, one can compare the thermodynamic functions of the 
probe with the variations of these same functions for the background when Nf is increased 
by one unit. For consistency, these two quantities should be equal. Actually, within the 
probe approximation one should assume that Nf is large. Then, the variation induced in the 
background when Nf — > Nf + 1 should be computed by a Taylor expansion in which only 
the first term is kept. We will verify that this compatibility condition is indeed satisfied in 
our case, which is a highly non-trivial test because the dependence of the background on 
Nf is non-linear. After passing successfully this test, we are ready to study systematically 
the thermodynamics of the probe brane. In general, the main objective is to determine the 
dependence of the different observables on the number of flavors of the background, as well 
as the departure from conformality induced on the system by the probe. 

The plan of the rest of this paper is the following. In Section [2] we will present our flavored 
black hole background and compute some of its thermodynamic functions. In Section |3] wc 
will analyze the flavor brane embeddings at zero temperature and extract some useful infor- 
mation which will be needed in the black hole case. In Section H] we will study the action 
of the probe in the non-zero temperature geometry and we will check that the compatibility 
condition mentioned above is satisfied. In Section [5] we shall study in detail the two types 
of embeddings, Minkowski and black hole, and we shall analyze the first order phase transi- 
tion between them. Section |6] is devoted to the calculation of the different thermodynamic 
functions of the probe (free energy, internal energy, entropy, and normal speed of sound). 
Section [7] contains a summary of our results and a discussion. The paper is completed with 
several appendices, which contain some explicit calculations and details not included in the 
main text. 
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2 The flavored ABJM background 



In this section we will present the non-zero temperature version of the deformed ABJM 
background found in [16]. The ten-dimensional metric, in string frame, of this supergravity 
solution takes the form 

ds 2 = L 2 ds 2 BHi + ds\ , (2.1) 

where L is the radius of curvature and ds 2 BHi is the metric of a black hole in the four- 
dimensional Anti-de Sitter space, given by 

ds 2 BHi = ~r 2 h(r)dt 2 + + r 2 [(dx 1 ) 2 + (dx 2 ) 2 ] , (2.2) 

and ds\ is the metric of the compact internal six-dimensional manifold^] In (12. 2p the black- 
ening factor h(r) is given by 



r 



h(r) = 1 - -A , (2.3) 



r 



where the horizon radius is related to the temperature by 



2tt 



d 



fg^ r dr 



3r h 



zhFto =-rr- (2-4) 



r=r h 



4lT 



The internal metric ds\ in (12 .ip is a deformed version of the Fubini-Study metric of 
This deformation is due to the backreaction of the massless flavors, generated by the D6- 
branes, and can be simply stated when the manifold CP 3 is represented as an S 2 -bundle over 
S 4 , with the fibration constructed by using the self-dual SU{2) instanton on the four-sphere. 
Explicitly, this metric can be written as 



L 



. 2 

dS l ~ ~p 



qds 2 s , + (dx' + e ijk A j x k ) 2 , (2.5) 



where b and q are constant squashing factors, ds^ is the standard metric for the unit four- 
sphere, x l (i = 1,2,3) are Cartesian coordinates that parameterize the unit two-sphere 
{^2i{x 1 ) 2 = 1) and A % are the components of the non-Abelian one-form connection corre- 
sponding to the SU (2) instanton. 

The squashing factors q and b in (12.51) encode the effect of the massless flavors in the 
backreacted metric. Indeed, when q — b = 1 the metric (12. 5 j) is just the canonical Fubini- 
Study metric of a CP 3 manifold with radius 2L and (12. ip is the metric of the unflavored ABJM 
model at non-zero temperature. The parameter b represents the relative squashing of the 
CP 3 part of the metric with respect to the AdS± part due to the flavor, while q parameterizes 
an internal deformation which preserves the S 4 -§ 2 split of the twistor representation of 1 



1 Unlcss otherwise stated, we will use units for which a! = 1. 
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The explicit expression for the coefficients q and b of the smeared solution of [16] is given 
below. They depend on the number of colors N and flavors Nf, as well as on the 't Hooft 
coupling A = N/k, through the combination 



e = 



3Nj 
Ak 



3 Nf_ 
A N 



A 



(2.6) 



where the factor 3/4 is introduced for convenience. The AdS radius L can be also expressed 
in terms of A and the deformation parameter (I2.6P (see eqs. (12.291) and (12.301) ). 

The type IIA supergravity solution found in [16] contains, in addition to the metric (12.11) . a 
constant dilaton and RR two- and four- forms F 2 and _F 4 . In order to specify the form of the 
latter, let us introduce a specific system of coordinates to represent the metric (I2.5|) . First of 
all, let u) % (i = 1, 2, 3) be the SU{2) left-invariant one-forms which satisfy duj 1 = ~ u 3 ' Au k . 
Together with a new coordinate £, the u ll s can be used to parameterize the metric of a four- 
sphere S 4 as 



2^2 



d e + 1 ((c 1 ) 2 + (c 2 ) 2 + (c 3 ) 2 ) 



(2.7) 



where < £ < oo is a non-compact coordinate. The SU(2) instanton one-forms A % can be 
written in these coordinates as 

e 



A' 



to 



(2.8) 



Let us next parameterize the x l coordinates of the § 2 by two angles 9 and (p (0 < 9 < it, 
< < 2tt), namely 



x 



sin 9 cos p , 
Then, one can easily prove that 



x 2 = sin 9 sin tp 



(dx l + e l3k A j A k ) 2 = (E 1 ) 2 + (E 



x ' 



cos 9 



2\2 



where E l and E 2 are the following one-forms: 



£2 

E 1 = d9 H — (sin ip tu 1 — cos ip ui 2 ) 

1 ~~\~ c 



(2.9) 



(2.10) 



(2.11; 



E 2 = sin 9[ dp 



e 



00 



e 



i+e j i + e 

Using these results we can represent the ten-dimensional metric (12.11) as 

L 2 



cos9 (cos^w 1 + sinpco 2 ^) . (2.12) 



ds 2 = L 2 ds 2 BH4 + . 



qdsl + (E^) 2 + (E 2 ) 2 



(2.13) 
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We shall next consider a rotated version of the forms of by the two angles 9 and ip. Accord- 
ingly, we define three new one-forms S l (i = 1, 2, 3): 



S* 1 = sin cp co 1 — cos p u 2 

S 2 = sin 9 to 3 — cos 9 (cos p tu 1 + sin </? w 2 ) 

S* 3 = — cos 9 u 3 — sin 9 (cos + sin p ui 2 ) . 



(2.14) 



In terms of the forms defined in (12.141) the line element of the four sphere is obtained by 
substituting u % —> S l in (12 .7p . Let us next define the one- forms S^ and S\ 



S' 



i + e i + e 

in terms of which the metric of the four-sphere is 

dsl = (^) 2 + 5>*) 2 



1,2,3) 



(2.15) 



(2.16) 



With these definitions, the ansatz for F 2 for the flavored background written in eq. (5.6) of 
ref. [TB] is 



k 



E 1 A E 2 - r] A S 3 + S 1 A S 2 ) 



(2.17) 



where rj is a constant squashing parameter between the § 4 and S 2 components of (I2.17p . In 
the unflavored ABJM solution of [3] the F 2 is given by (12.171) with rj — 1. For a general value 
of ?7 the two-form F 2 is not closed. Indeed, one can easily verify that 



dF 2 



2vr Q 



where Q is the following three-form 

n = ^- (1-r])^ 1 A(S^ AS 2 - S 1 AS 3 ) + E 2 A(S^ AS 1 + S 2 AS 3 ) 



(2.18) 



(2.19) 



Thus, when 77 7^ 1 the Bianchi identity for F2 is violated. This violation is due to the presence 
of a delocalized set of D6-branes, whose Wess-Zumino action can be written as 
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wz 



T, 



D 6 



C 7 AQ 



(2.20) 



Mi 



where C7 is the RR seven-form potential and Q is a charge distribution three-form. Clearly, 
the term (I2.20p induces a source for C7, which modifies the Maxwell equation of F$ = dCj. 
Taking into account that F 2 = *F S , one easily concludes that the equation of motion for 
Cj just takes the form of the modified Bianchi identity (I2.18p . Thus, one identifies the 
three-form Q written in (1 2 . 1 9 [) with the one parametrizing the distribution of the smeared 



set of D6-branes. Actually, from this identification one can relate the constant 77 to the total 
number of flavors Nf. Indeed, one gets [16] the simple equation: 

77 = 1 + ^, »7e[l,oo). (2.21) 

It is obvious from (12.211) that rj is simply related to the deformation parameter introduced 
in (J2SD, 

77 = 1 + e . (2.22) 

In the solution of [16] the squashing parameters b and 77 are related by a quadratic equation, 
which is obtained by requiring that the background preserves J\f = 1 supersymmetry at zero 
temperature. This quadratic equation is 

g 2 - 3(1 + 7])q +5r/ = . (2.23) 

By solving this equation for q and using (12.211) one can obtain q as a function of the defor- 
mation parameter e, 

d = 3 + - 2^1 + 6 + ^e 2 . (2.24) 



Moreover, the solution of the BPS system of [16] allows to relate the parameter b to the 
squashing factors q and 77: 

b = , q{v + q) r . (2.25) 
2(q + 77(7-77) 

From this equation we get the explicit expression of b in terms of the deformation parameter 
e: 

4+fe- Jl + e + i^ 

b = — . (2.26) 

3 + 2e v ' 

By construction 77 = q = b = 1 when Nf = 0, whereas in the flavored solutions these 

coefficients are greater than one. In order to have a better idea of the behavior of q and b it 

is quite useful to expand them in powers of Nf/k. We get 

, 3 N f 45 ,N f .2 , , 3 N f 63 /iV f \ 2 

q = l + -- L + •••, 6 = 1+ +•••• 2.27 

y 8 k 256 ^ fc ^ ' 16 k 512 V k J K J 

Notice, however, that q and b reach a finite limiting value when the deformation parameter 
is very large. Indeed, one can check from (12. 24ft and (12. 26ft that 

5 ,5 N f , 

q^-, b -)• - , as X^ 00, (2 ' 28) 

To fix completely the metric (12.11) we need to know the value of the AdS radius L. In the 
unflavored case L 2 is proportional to the square root of the 't Hooft coupling A. This value 
gets deformed by the backreaction of the flavors. Actually, we have 



L 2 = it V2X a , (2.29) 
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where a is defined as the following function of the deformation parameter: 



a = x j ? 2 ~ q , b> = \ ^ + g) 2 ( 2 -g) 5 . (2.30) 



It was shown in [16] that cr characterizes the corrections of the static quark-antiquark po- 
tential due to the screening produced by the flavors. In Fig. [1] we depict q, b, a, and £ (12.361) 
as functions of the deformation parameter e. 




Figure 1: Representation of the squashing factors q and 6, the screening function cr, and the 
volume function £ (I2.36P for the background, in terms of the deformation parameter e. 



The solution is completed by a constant dilaton given by 

-4, = b r] + q k 
4 2-q L ' 

and a RR four-form F 4 , whose expression is 



(2.3i; 



Fi = 3k L2 ^ (2 32) 

where Qbh 4 is the volume- form of the four- dimensional black hole (12. 2p . The regime of 
validity of the type IIA supergravity description can be obtained by requiring that L >> 1 
and e^<l. For the flavored ABJM background at zero temperature these two conditions 
were worked out in detail in [16] and will not be discussed further here. 

In the zero temperature case this background was found in [TB] by solving the system of 
first order BPS equations required to preserve M = 1 supersymmetry. Then, one can verify 
that the solution satisfies the second order equations of type IIA supergravity with sources 
(see appendix D of [IS])- In the black hole case one can easily check that these equations of 
motion are still satisfied after the introduction of the blackening factor h{r) in the metric. 



10 



2.1 Thermodynamics of the background 



Let us now find the values of the different thermodynamic functions for the flavored black 
hole presented above. We begin by computing the entropy density Sback, which is given by0 

2tt As 

Sback — 9 ] (Z.OO) 

where A 8 is the volume at the horizon r = r/, of the eight-dimensional part of the space 
obtained by setting r,t = constant in the ten-dimensional geometry and Vi is the infinite 
volume of the 2d space directions x l . The volume A$ has to be computed with the Einstein 
frame metric, which in our case is obtained by changing L by e - ^' 4 L in fl2.ll) and (j2.5p . 
After a simple calculation one can check that Ag/V^ is given by 

As _ 6*K_ q L e 3 

V 2 ~ 3 ^ r h- V-U) 

We can now use the values of the different factors appearing on the right-hand side of (I2.34p 
to obtain the value of the entropy density in terms of gauge theory quantities. Taking into 
account that, in our units, 2/cf = (27r) 2 , we get 

1 /4vr\ 2 N 2 _ fN f 



— -UtJ Wx t[ ^ ]T • (2 ' 35) 



£|'^ s ' '*'" + ') 4 ,. (2.36) 



where 

\ k J 16 ^/2^(q + 775 — 77)2 

The quadratic dependence of the entropy with the temperature is a reflection of the confor- 
mality of the system which, in our solution, is not affected by the massless flavors. Notice 

3 

that Sb ac k displays the characteristic TVs behavior of the effective number of degrees of free- 
dom of the ABJM theory in the 't Hooft limit. The correction to this behavior introduced 
by the flavors is parameterized by the function £, which was introduced in [16] and shown 
to be very close to the function obtained by using the localization technique. The function 
£ determines how the volume of the internal manifold (and, hence, the area of the horizon) 
changes due to the addition of flavor. 

The internal energy density can be obtained from the ADM energy, 



E 



ADM 



— vT^ttf / VdetGs (K T - K ) . (2.37) 

^10 JM t , roo 



In (I2.37P Gs is the Einstein frame metric of the t,r = constant hypersurface. The integral is 
taken over this hypersurface for a large value r = r ro of the radial coordinate. The symbols 



2 We use the same conventions as the first paper in [20] . 
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Kt and K denote the extrinsic curvatures of the eight-dimensional subspace within the 
nine- dimensional (constant time) space, at finite and zero temperature, respectively. For an 
arbitrary hypersurface K is given by 



K = -== d^^tG~ 9 n^ , (2.38) 

with n M being a normalized vector perpendicular to the surface. For a constant r hypersur- 
face, 

n" = -^=*?, (2.39) 
and one can show that K for our background becomes 

K = -— . (2.40) 

By using these results it is easy to find the value of the integrand in (12.37p . 

vlG^\/detG 8 (K T - K ) = -e^L 2 ^detg 6 r\ , (2.41) 



where g§ is the internal metric ( 12. ip . It is now immediate to obtain the internal energy 
density of the flavored black hole, 

E ADM 2 [Aix\ 2 N 2 j_ (N f 



Again, the dependence on the temperature is just the one expected for a conformal system 
and the flavor dependence is determined by the function £. Moreover, the free energy density 
Fback can be obtained from the thermodynamic relation F^ack = E^ack — T Sb ac k, yielding, 



back 



1 /4ttV N' 2 /N f \ 

U-r ) T ■ ( 2 - 43 ) 



9 V 3 J V2\ V k 



As a consistency check we notice that Sback = —dFb ac k/dT, as it should. It is also worth 
pointing out that the free energy density F back can be computed directly from the regularized 
Euclidean action (see the first paper in [20] for a similar calculation for the D3-D7 black 
hole). The regularization is performed by subtracting the action at zero temperature with 
the Euclidean time suitably rescaled. Furthermore, in the action one must include the 
standard Gibbons- Hawking surface term. The final result of this calculation, which will not 
be detailed here, is just the same as in (I2.43p . 
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3 D6-brane embeddings at zero temperature 



One key objective of this paper is to study the properties of flavor brane probes embedded in 
the flavored black hole background described in Section [2] Before dealing with this problem 
in full generality, let us analyze the case in which the temperature of the background is 
zero, which corresponds to taking the blackening factor h(r) equal to one in the formulas of 
Section [2j 

The kappa symmetric embeddings of the flavor D6-branes that preserve the supersymmetry 
of the zero temperature background were studied in [16] . As argued in [10] , these D6-branes 
should extend along the three Minkowski directions x^, the radial coordinate r, and wrap 
a three-dimensional submanifold of the compact internal space. For large values of the 
radial coordinate the metric of this three-dimensional submanifold should approach that of 
a (squashed) MP 3 = § 3 /Z 2 . In our § 4 — S 2 representation, it was shown in [TH] that this 
internal submanifold is obtained by extending the D6-brane along the § 4 base in such a 
way that the pullback of the one-forms u 1 and u 2 vanish. Accordingly, let us consider a 
configuration such that Cj 1 = Cj 2 = 0, where the hat denotes the pullback to the D6-brane 
worldvolume. Moreover, for the pullback of u 3 we just take u> 3 = dip, where ip is an angular 
coordinate. We will also assume that the brane is extended along the coordinate tp of the § 2 
fiber and that the other S> 2 coordinate 9 is a function of the radial coordinate r, 9 = 9{r). 
Therefore, we will choose the following set of worldvolume coordinates 

C = {x»,r,t,j>,<p) . (3.1) 

Then, the induced metric (at zero temperature) on the D6-brane worldvolume becomes 



ds 2 = -L 2 r 2 dt 2 + L 2 r 2 [ (dx 1 ) 2 + (dx 2 ) 2 } + L 2 
where ds\ is the following three-dimensional metric 

e 



1 1 /dB\* 
r 2 b 2 \ dr 



AL 2 

dr 2 + —ds 2 3 , (3.2) 



ds\ 



q 



+ 7 



(i+a 2 

If we redefine the angular coordinates as 



| ( 1 + ^ 2 )2 ^ 2 + £ sin 2 0( dp - Y^2 d $y 



tan 



a 



P 2 ' 



iff = cp 



then the 3d metric ds 2 becomes 



3 4 



(3.3) 



(3.4) 



(3.5) 



qda 2 + q sin 2 ad/3 2 + sin 2 (dip + cos ad/3)' 

where 9 is assumed to be a function of r. The range of the angular coordinates in (13. 5p is, 

0<a<7r, 0</3<2tt, < ^ < 2vr . (3.6) 
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Notice that, in these coordinates, the massless configurations whose backreaction is in- 
cluded in the background of Section [2], correspond to embeddings with 9(r) being constant 
and equal to tt/2. In order to simplify the study of all possible embeddings that satisfy the 
equations of motion of the probe, it is convenient to choose an isotropic system of coordi- 
nates. To find these coordinates, let us consider the (r,9) part of the induced metric (I3.2p . 
which can written as, 



dr 2 



L 2 



d9 2 



dr 2 + d9 z 



(3.7) 



We want to find a new radial coordinate u such that the first term inside the brackets in 
(13. 7p becomes du 2 /u 2 and the whole right-hand side of (13. 7p is proportional to du 2 + u 2 d9 2 . 
Clearly, we must require 

b JL = *L (3.8) 
r u 

and thus (13. 7p becomes 



b 2 v? 



[du 2 + u 2 d9 2 ] 



Eq. (13 .8p can be immediately integrated, with the result 

u = r b . 



(3.9) 



(3.10) 



Notice that the change r — > u of the radial coordinate is only non-trivial in the flavored case 
with b ^ 1. In terms of this u variable, the ten-dimensional metric (12. ip (for h = 1), becomes 



ds 2 



2 2 1 du 2 



dsl , 



(3.11) 



where dsl * s the metric (12. 5p of the squashed 

Let us now introduce a system of Cartesian-like coordinates (p,R), defined as 

R = u cos 9 , p = u sin 9 . 

The inverse relation is 



R 



u 2 = R 2 + p 2 , tan6> 

and, since du 2 + u 2 d9 2 = dp 2 + dR 2 , the line element ( 13. 7p becomes 

L 2 



[dp 2 + dR 2 ] . 



(3.12) 
(3.13) 

(3.14) 



b 2 (p 2 + R 2 ) 

Let us now consider embeddings of the D6-brane in which R = R(p). Then, the induced 
metric takes the form 



ds 2 = L 2 [p 2 + R 2 Y dx\ 2 + 



L 2 l + R 



12 



b 2 p 2 + R 2 



dp 2 + 



-2 r 



+ 



b 2 



q da 2 + q sin 2 adfi 2 + 



ir 



R 2 



[dip + cos a dp)' 



(3.15) 
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with R' = dR/dp. The embeddings corresponding to massless flavors are the ones for which 
R = 0. In the general case, the determinant of the induced metric takes the form 



L 



7 



y/-detg 7 = — q sinap[p 2 + R 2 }^- 1 Vl + R' 2 . (3.16) 

In order to obtain the explicit form of the embeddings, let us now study the action of the 
probe brane. We begin by computing the DBI action, which is given by 



Sdbi = ~T D6 J d 7 Ce-f V-det 



97 , (3.17) 



where the tension of the D6-brane Tdq = l/(27r) 6 in our units. Let us use ( 13 . 1 6[) in (I3.17P 
and integrate over the angular coordinates a, (3, and ip. We define a Lagrangian density 
£-dbi as 

Sdbi = / d 3 xdpC DBI , (3.18) 



where 

Cdbi = -No p[p 2 + R 2 }^- 1 Vl + R 2 , (3.19) 
with Mq being the following constant 

A/o = — p — - T m e~* . (3.20) 
Next, let us compute the WZ term of the action, which becomes 

Swz = Tub / C 7 , (3-21) 



where C-j is the RR seven-form potential (F$ = dCj) and, as before, the hat denotes the 
pullback to the worldvolume. In this zero temperature case the RR seven-form potential CV 
is naturally given in terms of the calibration seven-form K, that characterizes the G-structure 
of the supersymmetric solution. Indeed, we can take CV as 

C 7 = e-+K. (3.22) 

The seven-form K. is naturally defined in terms of a fermion bilinear which, in turn, can 
be obtained from the projections satisfied by the Killing spinors of the background. This 
calculation was performed in [16] and here we will limit ourselves to recall this result. As 
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shown in [16], to represent /C it is useful to define the following basis of one-forms: 
e° = Lr dt , e 1 = Lr dx , e 2 = Lr dy , 



e 3 = -rfr, e 4 = ^^, 

r o 

e* = f v^ 1 " 4 , (i = 5,6,7), 

e J ' = ^^ 7 , (j = 8,9), (3.23) 

which are a frame basis for the zero-temperature version of the metric ( 12. ip . In terms of the 
forms (13.231) the form K. can be written as [T6j . 

JC = -e 012 A (e 3458 - e 3469 + e 3579 + e 3678 + e 4567 + e 4789 + e 5689 ) . (3.24) 

To evaluate the WZ action we need to compute the pullback of K to the worldvolume. 
Let us write the pullbacks of the frame one-forms (13.231) in the (p, R) coordinates. In this 
calculation it is convenient to use 

R-pR' J 1 RR' + p A 
dO = dp , dr = — dp . (3.25) 

P + R b [p 2 + i? 2 ] 1 "^ 



We find 



= L[p 2 + R 2 ]^dx^ , e 3 = \ R f + J dp , e 4 = \ y/q da , 

b p z + R 2 b 

e 5 = , e 6 = — r~ sin a — , = : dB , e 7 = sin a — . = : dB , 

b y/p 2 + R 2 b ^p 2 + i? 2 

e 8 = \ ^§dp, e 9 = f ^=(# + cos^) . (3.26) 

By inspecting these pullbacks one readily verifies that the only non-zero contributions to K, 
are 

K, = e 012 A ( e 3469 - e 4789 ) = ^ sinap [p 2 + R 2 }^ 1 d 3 x A dp A da A dB A dip . (3.27) 
Thus, after integrating over the angular variables, we can write 



Swz — I d xdpCwz , (3.28) 
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with the Lagrangian density 

C wz = M oP [p 2 + R 2 Y b ~ l . (3.29) 
Therefore, the total Lagrangian density is 



C = -Nop [p 2 + R 2 ] * 1 ( VI + R' 2 - 1 ) . (3.30) 

Clearly, R = constant is a solution of the equations of motion derived from £ (notice that the 
on-shell action for this configuration vanishes). This is just the kappa symmetric solution 
that preserves SUSY which was found in [16]Jf| Let us now study the form of a general 
solution in the UV region of large p. In this case one can approximate p 2 + R 2 ~ p 2 in 
(I3.30p and take R' small. At second order in R', we find that C can be approximately taken 
as 

£ ~ _^ p !-ii2'2 . (3 . 31) 
The equation of motion derived from this second-order Lagrangian is simply 

^(pf- 1 it") = 0, (3.32) 

and can be integrated trivially 

R ~ m + -4- , (3.33) 
P~ b 

In ( I3.33P m and c are constants, which should be related to the mass of the quarks and 
to the vacuum expectation value of the corresponding bilinear operator ijjif) (see below), 
respectively. The power of p of the subleading term in (13.331) should determine the conformal 
dimension of the bilinear operator. Indeed, let us consider a canonically normalized field 
in Ad S4 with conformal dimension A. The behavior of near the boundary of AdS^ is 



or A " 3 + ^, (3.34) 



(O) 



r 



where 0o (the boundary value of 0) is identified with the source of the dual gauge theory 
operator O and the coefficient (O) is identified with its VEV. In (I3.34p A is the dimension 
of O and r is the canonical coordinate of AdS^ (in terms of which the AdS^ metric takes the 
form r 2 dx\ 2 + dr 2 /r 2 ). It is clear that this canonical coordinate is just the one in f !2.2p . In 
the UV, r and p are related as r ~ p 1 ^, and therefore we can rewrite f !3.33|) in terms of r as 

R-m+-^ b . (3.35) 

3 Notice that, in the angular (r,8) parameterization of [TB], the R = constant solution reads 8(r) = 
arccos (t 1 ) 6 , where r§ = Rq- 
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In order to extract the dimension of the operator dual to the scalar R, let us rewrite (I3.34p 
in such a way that the asymptotic value of the right-hand side is constant, 



(O) 



r d -^~0o+ ^A=3 ■ (3-36) 

Clearly, by comparing (I3.36P and ( I3.35P we find that, in our flavored ABJM case, 2A — 3 = 
3 — 2b, which yields 

A = 3 - b , (3.37) 

in agreement with the value obtained in [IB] for the dimension of the bilinear operator ipip. 
Notice also that A m = 3 — A is the dimension of the source (the mass in our case). This 
dimension is just A m = b in the flavored ABJM case. Thus, the mass anomalous dimension 
is 

7m = A m - 1 = b - 1 . (3.38) 

It is evident from f !3.38|) that the anomalous dimension 7 m depends on the number of flavors 
Nf and, according to (j2.28p . it becomes maximum when Nf — > oo: 

lT x = \ • (3-39) 

As it was already mentioned, the asymptotic value m should be related to the quark mass 
m q . To find the precise relation let us consider a fundamental string extended from the 
origin to the point with R = Ro = m8itp = 0. The induced metric on the worldsheet of 
this string is 

ds\ = -L 2 R>dt 2 + — — , (3.40) 

whose determinant is 



detc/ 2 = — Rb- 1 . (3.41 ; 

b 



The Nambu-Goto action for this string is 



1 f rR=m £2 

Sng = -77- / dt / dRy/-detg 2 = -— I dtm^ . (3.42) 



2vr ./ J R=0 2vr 



The action per unit time should be identified with m q . Thus, by using f)2.29p we arrive at 

a i ( m q \/a' \ b 



m oc — =■ m b =^r- m oc ( — ■= — ) , (3.43) 

where A is the 't Hooft coupling and a is the function of Nf/k that has been defined in 
(I2.30p . We have included a factor of \fa' = l s to reinstate the correct dimensions. 

The constant c in (I3.35P should be related to the vacuum expectation value of the meson 
operator ijjip (the quark condensate). In order to find this relationship we should relate c to 
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the derivative of the action with respect to the mass parameter m. In principle, to perform 
this calculation we should holographically renormalize the action to ensure its finiteness 
[2nj[27]. It turns out, however, that the action corresponding to the Lagrangian density 
f)3.30p is convergent and, therefore, this renormalization is not needed. Indeed, by using the 
asymptotic behavior f !3.35|) we obtain for large p, 



C ~ p x ~t , (3.44) 

and, since the maximum value of 1 — | is — |, the integral over p is convergent as claimed. 
Notice that this convergent behavior is a consequence of the particular gauge for Cj chosen. 
Indeed, performing a gauge transformation of the type C 7 — > C 7 + dA e is equivalent to 
adding a boundary term to the action of the probe and to choose a particular renormalization 
scheme. In our gauge Cj is chosen to be the calibration form and, as a consequence, the action 
for a supersymmetric embedding R = constant vanishes. For a more general embedding the 
WZ term introduces a subtraction of the DBI term, which renders the total action finite. 

The probe configuration is obtained by solving the equation of motion derived from the 
Lagrangian density (13.30P for R(p). In this process we have to impose boundary conditions 
at some value of the p coordinate. The simplest thing is to take p = as this initial value of 
the coordinate and to integrate the system outwards. It is easy to verify from the limit of the 
differential equation at p = that the only possibility to have non-singular solutions is to 
take R(p = 0) = R and R'(p = 0) = as initial conditions. At the UV region of large p the 
function R(p) must behave as in (I3.33p . where the constants m and c are not independent 
since both should be determined by the IR value Ro of R(p). The on-shell action is obtained 
by evaluating the integral of L for these configurations. It can be considered as a function of 
the mass parameter m. The derivative of S with respect to m can be computed as follows: 



dS 



dm 



dC dR dC dR 



OR dm OR' dm 



^ ~ d P ^^ + ^7^ = dP IT > ( 3 - 45 ) 







dp 



dC dR 



dR' dm 



where we have integrated by parts and used the equations of motion of R(p). In (I3.45P 
we have already integrated over the Minkowski coordinates and we have assumed that this 
integration gives rise to a constant factor. The value of the right-hand side of (I3.45P can 
be obtained by evaluating the "momentum" density dC / dR' at the boundary values of the 
worldvolume. It is readily checked that, for regular embeddings, the IR contribution at p = 
is zero. To obtain the UV contribution at p = oo, let us use (13. 35ft : 

H cx p [p 2 + R* ] v- 1 ~ 6 —^- c + subleading , (3.46) 

where we include different factors coming from the constant Mq. Taking into account that 

dR 

— — = 1 + subleading , (3.47) 
dm 
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we get 

95 3-21 . a \ 

c < 3 - 48 > 

The quark condensate (ipip) is obtained by performing the derivative of the action with 
respect to the bare quark mass fi q . The latter can be obtained by taking o = b = 1 in the 
dressed mass m q . It is clear from (13.431) that fi q ~ m and, thus, 

/7( . dS 3-26 , > 

w '«~fK~— c ' (3 - 49) 

Therefore, c is indeed proportional to the quark condensate. It turns out, however, that the 
only regular solutions in this T = case are those for which R = constant = m, i.e., the 
kappa symmetric ones. They have c = and therefore the quark condensate vanishes in 
this case. Notice that the on-shell action for these solutions is zero, as expected on general 
grounds from their supersymmetric character (see (13.3 



4 Flavor brane probes at non-zero temperature 



In this section we come back to the analysis of brane probes in the general non-zero tem- 
perature background of Section [2J The main difference from the T = analysis of Section 
[3] is due to the presence of an event horizon in the metric. Thus we will have two types 
of embeddings: Minkowski and black hole. In the former type the brane probe does not 
reach the horizon, whereas in the latter case the brane ends on the horizon. In order to 
describe correctly the thermodynamics of these two types of configurations and of the phase 
transition connecting them one has to define carefully the action of the probe. It turns out 
that there is a subtlety which we shall address in this section. 

As in the T = case, we will consider D6-brane probes embedded in the internal CP 3 in 
such a way that the one-forms u 1 and Gj 2 vanish. We will take (x M , r, a, f3, ip) as worldvolume 
coordinates and describe the embedddings by a function 9 = 9{r). The induced metric takes 
the form 



ds 2 = -L 2 r 2 h{r)dt 2 + L 2 r 2 [{dx 1 f + {dx 2 ) 2 ] + 



L 2 



r 2 h(r) 



r 2 h(r) 



e 2 



dr 2 + 



b 2 



q da 2 + q sin 2 a df3 2 + sin 2 6 (dip + cos ad(3)' 



(4.1) 



where the dot represent the derivative with respect to r. The determinant of the incuded 
metric is 

V q 



V-det 
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b 3 



r 2 sin 6 sin a \ / 1 H — — 6 2 



(4.2) 



20 



After integrating over the internal space we get the following DBI action: 



S DBI = M r I d 3 xdrr 2 sin^/l + 6 2 , (4.3) 



with M r being the following constant 



& 3 

In terms of gauge theory quantities, we can write M r as 



K = ^ T„ e~* . (4.4) 



AV2n Vk V k J 



where the function ((Nf/k) contains all the dependence on Nf and is given by 



N f \ 1 V2^(77 + g)6 4 1 V^Qiv + Q^q 7 



2 \Zq(<i + - r i) 2 ^2 (q + vq~v) 



(4.6) 



Notice that £ = 1 for iV/ = and for an arbitrary number of flavors this function is related 
to the screening function a defined in (I2.30P by a simple equation 

* = £C (4-7) 

Let us now focus on the WZ action, which requires some extra consideration to eventually 
yield consistent thermodynamics. Recall that the WZ term of the probe action is propor- 
tional to the integral of the pullback of the RR seven- form potential CV (see fl3.2ip ). In the 
zero-temperature case analyzed in Section [3] we represented CV in terms of the calibration 
form /C (eq. (13.221) ). Actually, one can easily verify that introducing the blackening factor 
h(r) does not change the field strength F 8 = — * F 2 (the dependence on h(r) cancels when 
one computes the Hodge dual of F 2 ). We verified in Section [3] that this is a gauge choice that 
leads to an on-shell action of the probe which is finite at the UV. Since in this region the 
modification of the background due to the temperature vanishes asymptotically, it is clear 
that Cj for T ^ should also contain /C. Moreover, in the general case we should worry 
about the behavior at the horizon. Let us explore the possibility to improve the behavior of 
the worldvolume action at the horizon without modifying its regular character at the UV. 
In general, we will write CV as 

CV = e~* K + 5C 7 , (4.8) 

where 5Ci is a closed seven-form which must vanish in the SUSY (zero temperature) case. 
To determine the improving term 5Cj in (14. 8p . we first study the pullback of K, in the black 
hole case. The expression of /C is the one written in (13.241) . where the e a are the one- forms 
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defined in f!3.23j) (notice that they do not contain the blackening factor). Recall that the 
angular embedding of the D6-brane is characterized by the conditions Cj 1 = Gj 2 = 0, which 
imply that e 5 = 0. Along this submanifold, the pullbacks of the one-forms in ( I3.23P are 



Lrdx^ 



(^ = 0,1,2) 



dr 



L 



y/qda 




L 



— r — sin a sin^c^ 



-— Jq sin a cos 9 df3 
b 



dQ 



sin 9 (dtp + cos a df3) 



(4.9) 



where the hat over the forms denotes the restriction to the angular submanifold defined by 



the conditions Cj 1 
is given by 



0. Using these results we get immediately that the pullback of /C 



L q j. r r 3 



If' 



e~* d 3 x A 



— sin 9 cos 9 d9 + r 2 sin 2 9dr 
b 



A 3 



3 • 



with 3 3 being the following three-form: 



—3 



sin a da A d/3 A dijj . 



(4.10) 



(4.11) 



Let us now represent the improving term SCi in a way similar to the right-hand side of 

HDD, 



5C 7 



L 7 q 



-* d 3 x A 



L 1 (9)d9 + L 2 (r)dr 



A3, 



(4.12) 



with L\{9) and L 2 {r) being two functions to be determined. Notice that 5Cj is closed when 
L\ is only a function of 9 and L 2 only depends on r. The pullback of the total C 7 takes the 
form 



Vq 
b 3 



^ d 3 x A 



— sin^ cos^ + L 1 (9))d9 + ( r 2 sin 2 # + L 2 {r))dr 



A3, 



(4.13) 



As argued in [13] (see also [E]), a non-zero value of CV at the horizon introduces extra 
sources in the theory which change the boundary conditions of the fields and should be 
avoided. Accordingly, we impose the condition that the angular part of CV (i.e., the one that 
does not contain dr) vanishes at the horizon r = r^. This regularity condition determines 
uniquely the function Li(6 l ), 



sin# cos 9 . 



(4.14) 



Notice that, for this value of -£q(0), one can recast the d9 component of Cj in terms of the 
blackening factor, 

(4.15) 



3 3 

— sin# COS6 1 + Li(6 ) ) = —h(r) sin# COS6 1 . 
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It is important to point out that this term always vanishes at the bottom of the brane which 
is either at r = (for black hole embeddings) or at 9 = (for Minkowski embeddings). 
Thus, the pullback of CV to the submanifold with u 1 = Gj 1 = is 



C 7 



L 7 q 



d 3 x A 



— h(r) sin 9 cos 9d9 H ( r 2 sin 2 t 



A H 



3 • 



and the WZ term of the action is given by 



rh(r) 



cos99) + Af r \ d 3 xdrL 2 {r 



Let us now introduce a constant An, defined as 



J drL 2 {r) = r 3 h A 



(4.16) 



(4.17) 



(4.18) 



where the factor r\ has been introduced for convenience and the definite integral is over the 
whole range of the radial coordinate. Then, 



Swz = -N't J d 3 x dr r 2 sin 9 ^ sin 9 



rh(r) 



cos 9 9 



M r rl \ d 3 xA 



(4.19) 



Clearly, as the constant An does not depend on the embedding, it is a counterterm that 
represents a zero-point energy]^ The total action is given by 



S = -M r / d 3 xdrr 2 sin 9 



¥ b J J 

(4.20) 

Notice that the canonical momentum for the improved action f!4.20p vanishes at the horizon, 



dC 

d9 







(4.21) 



r=r h 



This means that the IR contribution to on-shell quantities like the one in (I3.45P will vanish 
for black hole embeddings that end on the horizon. This is related to the fact that, due to 
(I4.2ip . there is no momentum flow through the horizon and thus the latter is not a dynamical 
surface. This property will be important in what follows. 

Let us now fix the zero-point constant An in f)4.20p . With this purpose we will compute 
the free energy of the probe and compare this result with the free energy of the flavored 
background that was obtained in Section |2J In general, the free energy F is obtained from 
the Euclidean action Se by the relation F = T Se- In the calculation of Se one has to 
integrate over the Euclidean time r in the range < r < 1/T and over the non-compact 



4 Howcver, its contribution to the free energy is not a thermodynamic constant since it is multiplied by 
T 3 , as it is clear from the factor multiplying it in (|4.19[) (see below). 
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two-dimensional space. The latter gives rise to an (infinite) two-dimensional volume V%. In 
what follows we will divide all extensive thermodynamic quantities by V2 and we deal with 
densities. In particular, the free energy density (which we will continue to denote by F) 
corresponding to the probe action (I4.20p is 



F = M r I drr 2 sm.9 



r 2 h(r) 



9 2 - sinO - 



rh(r) 



cos 6 6 



Kr 3 h A c 



(4.22) 



In the next subsection we will determine the constant A by considering the case in which 
the probe brane remains very far from the horizon. This case corresponds to having quarks 
with very large mass which should decouple and therefore should not contribute to the free 
energy. As we will soon demonstrate, the condition F(m q — > 00) = will determine a simple 
value for A . 



4.1 Decoupling infinitely massive flavors 



To characterize the embeddings which correspond to flavors with infinite mass it is very 
convenient to work in a system with isotropic (Cartesian-like) coordinates. Let us proceed 
as in the zero temperature case and find a coordinate u such that the (r, 9) part of the metric 
is written as in (13. 9p . It is immediate to conclude that, in this black hole case, the differential 
equation for u(r) is 

(4.23) 



bdr 
r\fh 



du 
u 



which again can be integrated straightforwardly 



_3_ 



V 



f \ 3 

r h 



Notice that the horizon r = ry l corresponds to u — 1. The inverse relation is 



r 



U 2b -\- U 2b 



-u* f{u) 



(4.24) 



(4.25) 



where we defined a new function f(u), 

f{u) = 1 + u-f . (4.26) 
Let us next define a function f(u) as 

f(u) = 1 - u~t> . (4.27) 
One can verify that the blackening factor can be written in terms of / and / as follows 

\fh = L. (4.28) 
/ 
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Let us next write the ten-dimensional metric of the ABJM flavored black hole in terms of 
the isotropic coordinate u. We have 



L 2 r 2 



ds 2 



2l 



2 ~4 
3 



Cdt 2 + (dx 1 ) 2 + (dx 2 ) 2 

f 2 



+ 



b 2 u 2 



ds% , 



(4.29) 



where dsfl is the squashed CP 3 metric written in (12.51 ). 

Let us now define new coordinates R and p as in (I3.12p and parameterize the embedding of 
the probe by a function R = R(p). Following the same steps as above we can readily obtain 
the action of the probe and the corresponding free energy. The details of this calculation 
are given in Appendix |A] The total action for an arbitrary value of A is written in (lA.lip . 
By studying the p — > 00 limit of the equation of motion derived from (lA.llj) it can be easily 
proven that the function R(p) has the asymptotic behavior displayed in (I3.33P and therefore 
the solutions are characterized by two constants m and c, which are related to the quark 
mass and condensate, respectively. Moreover, from flA.llj) it is immediate to obtain the 
expression for the free energy density F. To write this result it is quite useful to define a 
new quantity Af as 

N =^ r h~ fti T ^ e ■ 

In terms of gauge theory quantities M has the following expression 



(4.30) 



Af 



2\/2 7r 



(4.31) 



where £ is the function of Nf/k defined in (14. 6p . Then, the free energy density for an 
embedding characterized by a function R(p) is 



F = Af 



dpp[p 2 +R 2 ] 2 " 1 ff Vl + R' 2 -1 



R 



p 2 + R 2 



(pR'-R) 



-AbAr 



, (4.32) 



where R' = dR/dp. This expression simplifies greatly when we take R = Rq 
this case we have 



FjR = R Q ) 
Af 



2R 2 Q 



dp 



P 



(P 2 



Rl) 2 



1 



(p 2 



R 2 )t b 



4b Ac 



constant. In 



(4.33) 



The integral on the right-hand side of (14.331) can be integrated straightforwardly. The result 
is 

F(R = Ro) 



1 - 



2b 



Rn 



46 A . 



(4.34) 



Af 2b + 3 

By looking at the equations of motion of the probe in the (R, p) variables it is easy to 
convince oneself that R = Rq = constant is a solution only in the case for which R — > 00, 
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which corresponds to the case for which the quark mass parameter m is very large. In this 
limit (14.341) becomes 

lim F ( R = R ^ = i _46A . (4.35) 

As argued above, infinitely massive flavors can be integrated out and therefore their contri- 
bution to the thermodynamic functions should vanish. Thus, on physical grounds one should 
choose A in such a way that the right-hand side of (I4.35P vanishes, which means that A is 
simply given by 

A = I . (4.36) 
4.2 A highly non-trivial test 

Let us show that the value of A written in (I4.36P is precisely the one required to satisfy a 
non-trivial compatibility condition between the free energy density of the probe and the one 
obtained from the flavored geometry. With this aim let us determine A again by considering 
the case of zero mass embeddings (for which 9 is constant and given by 9 = 7r/2). One can 
readily verify that this configuration solves the equations of motion derived from the action 
f)4.20p and that the only contribution to the free energy (14.221) is precisely given by the 
zero-point term. Thus, in this case we have 

/Aix\ 3 

F « -Mr r 3 h A = - ( T ) Mr A T 3 . (4.37) 

At this point it is interesting to remember that our background includes the backreaction of 
Nf massless flavor branes. Actually, the free energy (I2.43P contains the effects of Nf flavor 
branes at non-linear order in Nf. In the limit of small mass the free energy of the probe 
should match the variation, at linear order, of the free energy of the backreacted background 
( I2.43P when one flavor is added. Let us compute this variation at linear order by expanding 
the function £ (jk) defined in (I2.36P in a Taylor series and keeping only the first order, 

where the prime denotes derivative of £ with respect to Nf/k. Therefore, the variation of 
the free energy (12.431) of the background (at linearized level) is 



By equating AF back with the right-hand side of ( I4.37P we find the following value of A 



A f 



Ni 1 (Nf 



12^71 fc3 M r V k 



CU) . (4.40) 
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To simplify this expression of Ao, let us rewrite J\f r as in (14. 5p . Then, we can readily check 
that all the dependence on N and k drops and the expression for the zero-point constant Ao 
is greatly simplified. We arrive at 

Ao = J£ • (4.41) 

Remarkably, by computing explicitly the derivative with respect to the deformation param- 
eter Nf/k = e one can find a simple expression of £' in terms of q, t], and b for arbitrary 
values of the deformation parameter. This expression is 

!, for^->0, 



2 



512 y 2 A/ JV) ' iUI fc ^ °° 



where the limiting cases match with eqs. (7.9) and (7.11) of [16], respectively. Amazingly, 
this value of £' is simply related to the function £ that encodes the flavor dependence of the 
prefactor of the probe free energy. Actually, by comparing the right-hand sides of (I4.42p and 
( 14. 6 p one readily concludes that 

«'=|C, (4.43) 

which, after taking (14.41 p into account, means that A is just given by ( I4.36p . in remarkable 
agreement with our calculation in the opposite m — > oo limit. 

The result just found implies that the first variation of the free energy of the flavored back 
hole can be written as 

AF back = ~^r\. (4.44) 

Obviously, in terms of A/", the first flavor variation of the free energy of the background takes 
the form 

AF back = -M. (4.45) 
It follows that the limiting value of the free energy for massless embeddings is 

F w -N , asm^O. (4.46) 

It is interesting to formulate the matching between the action of the probe and background 
in terms of the entropy density. In the m — > limit the entropy density of the probe is just: 

s = ~%^-!fi (m->0), (4.47) 

which, after using ( I4.3ip and (I4.43p . can be written as: 
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Let us now calculate the total entropy of the system, i.e., the sum of (14.481) and the back- 
ground entropy (I2.35|) . By linearizing the function £ as in (|4.38|) . we can write 

1 /4tt\ 2 N 2 ^ (N f + 1\ - , . , > 

stotai = s back + s « - — -=d — 7 I T , (m -> 0) , (4.49) 



2A 

which means that Stotai is equal to the entropy of the flavored black hole in which Nf is 
increased by one unit. Therefore, the effect of adding a probe brane with m — > is equivalent 
to the increase of the area of the horizon which is produced in the geometry when Nf — > Nf+1 
and, thus, the effect of the probe in this limit is very nicely encoded in the geometry of the 
backreacted background. 

Notice that the dependence on Nf of the entropy of the background is determined by 
the volume of the squashed CP 3 manifold, while that of the massless probe is related to 
the volume of the squashed 1RP 3 cycle that it wraps. Thus, the compatibility condition just 
checked means that the volume of the cycle is simply related to the derivative with respect 
to Nf of the total volume of the internal manifold. Given the fact that these volumes depend 
non- linearly on Nf, this is a remarkable property of the background which we regard as a 
highly non-trivial test of the consistency of our flavored geometry. 



4.3 Summary of the RR potential and action 



To finish this section let us summarize the result of the previous discussion. We have found 
that the RR seven-form potential which satisfies the requirements imposed by the holographic 
renormalization and regularity at the horizon of the flavor brane must have the form: 

K 



c 7 



X + 



8tt 2 T D6 



d 3 x A 



l(r) dr — 4 sin 9 cos 9 d9 



A S 



3 , 



(4.50) 



where K. is the calibration form (I3.24p . M is written in (I4.30p . S3 is the three-form (14. lip 
and l(r) is a function whose integral over r must be one in order to decouple the infinitely 
massive flavors. If the embedding of the brane is parameterized by a function 9(r), the total 
action of the probe is given by 



s = K 



d 6 x 



1 it I drr 2 sin 9 



r 2 h(r) 

~9 r ~ 



9 2 - sin 9 



rh(r) 



while in terms of the (R, p) variables becomes 



S 



Af / d 3 x 



dpp[p 2 + R 2 ] 2b 1 ff Vl + R' 2 - 1 



_R 

7 



R 2 



cos 9 9 



(pR' - R) 



(4.51) 



(4.52) 

Once the action is completely fixed it is rather straightforward to study the different 
solutions of the equations of motion and their corresponding thermodynamical properties. 
This analysis will be carried out in the next sections. 
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5 Minkowski and black hole embeddings 



The action (I4.52p is certainly more complicated than its zero temperature counterpart (I3.30p . 
However, in the UV region of large p the equation that determines R(p) is still given by 
(I3.32p and therefore the embedding function R(p) behaves asymptotically as in (I3.33p . The 
constants m and c are related, respectively, to the quark mass m q and to the quark condensate 
{O m ). The detailed relation between m and m q is worked out in Appendix [Dj and is given 
by 

2^7T /— m 1 , 

m q = — - V2A a Tm* , (5.1) 

where A = N/k is the 't Hooft coupling and a is the screening function defined in f)2.30p . 
Notice that, according to fl5.ll) . taking m — > (m — > oo) for fixed m q is equivalent to sending 
T — > oo (T — > 0). Moreover, following the same steps as in the zero temperature case, we 
can relate the constant c to the quark condensate (O m ). Indeed, it is proved in Appendix [D] 
that this relation is 




Figure 2: Different embeddings in the (R, p) plane for the unflavored background (left) with 
m — 1.8 and for the flavored one with e = 10 (right) with m = 2.5. 

At low temperature (or large mass parameter m) the probe brane closes off outside the 
horizon and one has a Minkowski embedding. In this case the brane reaches the point p = 
(or 6 = 0) where the coordinate R takes the value R(p = 0) = Rq. One can readily check 
that the only solutions of the equation of motion derived from the Lagrangian flA.lip which 
are non-singular at the endpoint p = are those such that R' = 0. By imposing these two 
initial conditions at p = one can integrate numerically the equation of motion and find 
the function R(p). Some of these solutions for different values of Ro are shown in Fig. |2j In 
general the value of Rq determines the the asymptotic constants m and c and, by eliminating 
Ro, one can determine c = c(m). For general values of Rq this relation can only be found 
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numerically (see Fig. |3]). However, for large Rq (or, equivalently, large m or small T) one 
can establish an approximate relationship. Indeed, it is shown in Appendix iBl that 



m ~ R + 



3 + 26 



2b 



3-26 



Rn 



(Ro,m large) 



(5.3) 



where ip(x) = T'(x)/T(x) is the digamma function. Moreover, in this low T regime one can 
also obtain the function c(m) for large m, which is given by 



66 



46 2 - 9 m i+f ' 



m ^> 1 



(5.4) 




Figure 3: On the left we plot c versus m. The solid curve corresponds to the unflavored 
background while the dashed curve is for e = 10. In both curves the black color stands 
for black hole embeddings, while the blue for Minkowski. On the right we present a zoom 
showing the spiraling behavior near the phase transition point. 

When the temperature is large enough the probe brane ends at the horizon and we have a 
black hole embedding. In this case it is more convenient to use the isotropic coordinate u as 
the holographic coordinate and to represent the profile of the brane in terms of the function 
x(u), defined as 

x(u) = cos 6{u) . (5.5) 

The action in these variables has been obtained in Appendix |A] (see eq. flA.20j) ). The 
corresponding equation of motion for x{ u ) is 



0„ 



fx + 



V 7 ! - x 2 + u 2 x 2 



M- 71 = 2 . + 2 fx 

y/l - X +U 2 X 2 



, ( 5 - 6 ) 

where now the dot denotes differentiation with respect to u. From 05. 6 p we can infer the 
asymptotic behavior of x( u — > °°) : 

m c , . 

X — — H 3~ 7 + " ' ? 5.7 
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where m and c are the same constants as in (13.351) . Using the fact that f(u = 1) = one can 
immediately show that the solutions of (15.61) which are non-singular at the horizon u — 1 are 
those which satisfy the conditions 

X (u = 1) = X h , X(u = 1) = . (5.8) 

Some of the numerical solutions of (15. 6p with the initial conditions (15. 8 p are shown in 
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Figure 4: Phase transition temperature T c versus e. 

Fig. [2J In (I5.8p Xh is an IR constant which determines the UV constants m and c. As in 
the Minkowski embeddings, by eliminating Xh one gets c = c(m), a relation which can only 
be obtained for all values of m numerically. These results are plotted in Fig. [3j For high 
temperature (or low mass) x remains small for all values of u and one can linearize (15. 6p . 
which then can be solved analytically. This analysis is performed in detail in Appendix [C| 
where it is shown that, in this limit, xh is linearly related to m by 

r(i - ±) 

Xh « Vtt w \ _ jr m , m < 1 . (5.9) 

Notice that the coefficient multiplying m contains the dependence on the number of flavors. 
Similarly, one can find the function c = c(m) for small m, which is given by the following 
analytic equation 

c«- r(|)r(f-|) m ' m<<1 ' (5 - 10) 

which implies that c vanishes linearly as m — > with a slope that depends on the deformation 
parameter Nf/k. 

The temperature T c of the first order phase transition grows with the number of flavors 
as shown in Fig. HI This temperature is determined as the point where the curves of the free 
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energies of the black hole and Minkowski embeddings intercept each other. It is important 
to point out that the value of 1/m where the Minkowski-black hole transition takes place 
does not change much with Nf. Indeed, it (monotonically) decreases from being « 0.544 at 
Nf = down to 0.400 as Nf — > oo. However, from (15. ip we have — y/\ oc m~b a -1 , with 
a proportionality constant which does not depend on the number of flavors. This means 
that the flavor dependence of T c is dominated by the function er" 1 which, for large Nf, grows 
with the deformation parameter as This is precisely the behavior displayed in Fig. HJ 

The black hole and Minkowski embeddings are separated by a critical solution in which the 
brane probe just touches the horizon. This critical solution occurs for certain values m = m* 
and c = c* of the mass and condensate parameters. The detailed analysis of these critical 
embeddings is performed in Appendix [Fj where it is shown that they can be approximately 
represented near the horizon as R(p) ~ 1 + p. The solutions near the critical embedding 
display a discrete self-similarity behavior, as it corresponds to a first order phase transition. 
Indeed, as shown in Appendix [Fj the mass and condensate parameters exhibit an oscillatory 
behavior around their critical values and, as a consequence, the quark condensate is not a 
single- valued function of the mass. This last fact is clearly visible in the plots of Fig. [3j 



6 Brane thermodynamics 



In this section we address the main objective of this paper, the calculation of the different 
thermodynamic functions of the brane probe. The first of these quantities is the free energy 
density F, which can be obtained as in (I4.22p from the Euclidean on-shell action of the probe. 
Actually, the expression of F can be easily related to the integrals of the Lagrangian density 
C of ( lA.lip and (1A.20|) . Indeed, let V 3 be the value of the volume of three-dimensional 
Minkowski space and let us represent the on-shell Minkowski action of the brane in terms of 
a function Qim) by 

i con— shell 

= 1 - Q{m) . (6.1) 



V 3 Af 

Then, the free energy density F is given by 



A7 = 0{m) 



(6.2) 



The explicit expression for the function Q(m) can be straightforwardly obtained from the 
results of Section [5] For Minkowski embeddings parameterized by a function R(p), it is given 
by 



g(m) 



dpp[p 2 + R 2 ] 2b V/ Vl + R' 2 - 1 



R 



R 2 



(pR'-R) , (6.3) 
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while for black hole embeddings it is more convenient to represent Q(m) as 



duffu 



r-l 



u 2 x 2 



f 



(6.4) 



In (I6.3P and (16. 4 j) it is understood that -R(p) and x( M ) are the result of integrating the 
equations of motion with the regular boundary conditions at the IR which correspond to the 
UV parameter m. In Fig. |5] we plot F as a function of m _1 for both types of embeddings. 
Notice that the curves for Minkowski and black hole embeddings cross and show the typical 
"swallow tail" form, which is characteristic of first order phase transitions. It is important 
to point out that the improvement term (I4.14p that regularizes the behavior of CV at the 
horizon is essential in obtaining this behavior. 




F 
N 




Figure 5: We plot the free energy F/M versus 1/m for black hole (black curves) and 
Minkowski (blue curves) embeddings. The solid (dashed) curves are for the unflavored 
(flavored with e = 10) background. On the right an amplification of the phase transition 
region is shown. 

Let us now compute the entropy density s. We start from the definition of s as a derivative 
of the free energy F with respect to the temperature, which we organize as follows: 

8F . d (F\ F dj\f 

s = -&T = - M &T\Af) ~IfW (6 ' 5) 

Taking into account that Af ~ T 3 , we have 

dJ\f 3 



dT T 

and therefore we can write (16. 5p as 



Af , (6.6) 
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Let us now use the fact that for fixed quark mass m q the parameter m behaves as m oc T~ b 
(see (15. ip ) and thus |^ = —6^- Using the chain rule, the derivative appearing on the 
second term in (16.71) becomes 

d / F\ m d f F 



&TWJ ~ - h T^\M)- (6 - 8) 

The derivative with respect to m appearing on the right-hand side of (I6.8P has been computed 
in Appendix [D] (eq. flD.17j) ). By using this result, we can write 

Plugging the value of this derivative in (16. 7J) , we arrive at the following expression for the 
entropy s: 

s 3 F m , , 

_ = _____ (3-26) c. 6.10 

The first term on the right-hand side of (16.101) is the one expected in a system with conformal 
invariance in three dimensions for which F oc T 3 . The term in (I6.10p containing m and c 
represent the deviation from this conformal behavior due to the massive quarks introduced 
by the probe. Notice that it depends on the number Nf of massless quarks of the background. 
By using (16. 2p we can write s in terms of the function G(m), 

T-j^ = -3g(m) + 3 - (3 -26) cm . (6.11) 

In Fig. [H] we have plotted the numerical values of the entropy as a function of m _1 . We 
notice that s is always positive. As with the free energy, the regularization of CV at r = is 
essential to avoid having a pathological thermodynamic behavior for which s < for some 
values of m. 

We can also compute the internal energy E by means of the thermodynamic relation 
E = F + T s. Indeed, from (16.101) we get 

E = —2F - Af (3 - 2b) cm . (6.12) 

In terms of G(m), this expression can be rewritten as 

= -2G{m) + 2 - (3 -26) cm . (6.13) 

In Fig. [7] we plot E for different values of m~ l . 

We can now use the previous expressions and the numerical results to find the limiting 
values for the free energy, entropy and internal energy when m is small (or T is very large). 
Indeed, since Q(m) — > as m — > 0, it follows that 

F s E 

lim —. = -1 , lim T —. = 3 , lim —. = 2 , (6.14) 

m-)-0 AA m-tO AT m->-0 J\f 
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Figure 6: We plot the entropy Ts/M versus 1/m for black hole (black curves) and Minkowski 
(blue curves) embeddings. The solid (dashed) curves are for the unflavored (flavored with 
e = 10) background. 

which are just the values expected in this conformal limit. Moreover, in the opposite regime 
m — > oo (or T — > 0) one has Q{m) — > 1 and cm ~ m~ 3 ^ b — > 0. Thus, 

F s E 

lim — : = lim T —. — lim — : = . (6.15) 

TV m— »oo J\f m— >oo TV 



m— >oo 



In the next two subsections we will refine the limits (I6.14p and ( 16 . 1 5[) by using the results of 
Appendices [B] and O 

The heat capacity density c v of the probe is defined as 

OE 

cv = ■ (6-16) 

By computing the derivative of E as given in f)6.13p . one arrives at the following expression 
of c„: 



<9(log c) 



Tjj = 2Tjj - (3 - 2b) 3 - b - b Q ^_, cm. (6.17) 



d(\og m) 



We have checked numerically that c v is positive for all values of m and has a finite jump 
discontinuity at the phase transition point. 

6.1 Low temperature functions 

Let us now evaluate F, s, and E when T — > (or m — > oo). The on-shell action of the probe 
in this limit has been computed in Appendix [Bl From this result we find that Q (m) behaves 
as 

Gtrn) * I - . (6.18) 

ZO -f- O TTji, 
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Figure 7: We plot the internal energy E/Af versus 1/m for black hole (black curves) and 
Minkowski (blue curves) embeddings. The solid (dashed) curves are for the unflavored 
(flavored with e = 10) background. 



Moreover, the approximate value of c(m) in this T — > regime has been written in (15 .4p . 
By using this result in (I6.2p . (16. lip , and (I6.13p . we get 

F 2b fTy T s 12b /Ty E 10b /T\3 



Af 3 + 2b\MJ ' Af 2b + 3\MJ ' Af 2b + 3\MJ ' 

(6.19) 

where M is the constant defined in (ID. 31) . Taking into account that M ~ T 3 , we find that F 
vanishes as T 6 when T — >■ with a coefficient which depends on the number of flavors. As 
a check of (16.191) . one can immediately verify that the coefficients of F and s in (I6.19P are 
such that the thermodynamic relation s = —dF/dT is indeed satisfied. Furthermore, one 
can verify from (I6.19P (or directly from the general expression (I6.17P ) that the specific heat 
c v vanishes at low temperatures as T 5 , 

T$~^&)'. (6.20) 



Af 26 + 3 VM 
6.2 High temperature functions 

It follows from the results of Appendix O that Q(m) vanishes, when m — > 0, as 

. . 3 2b . 
y(m) ~ — — cm . (6-21) 
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Moreover, the condensate parameter c for small m can be estimated as in (]5.10[) . From these 
results we can show that F can be approximated as 



F 

Jr 



-i + 



r(i-t) 
r(l-l) 



My) (r) 



(6.22) 



from which it follows that the deviation of F from its conformal value decays as T~ 2h when 
T — > oo. Notice that, in this case, both the power of the temperature and the coefficient of 
this non-conformal contribution depend on the number of flavors. 

By combining (I6.1ip and (16.211) we can approximate the entropy in this limit as 



T 



3 - (3-26)<7(m) 



(6.23) 



which can also be written as 



T 



3(3 - 26) 



Af 



r(i-f) 
r(i-i) 



tan( y ) ( T ) 



(6.24) 



As a check of (I6.24p one can verify that s = —OF/ dT. Moreover, from ( I6.22p and (I6.24p we 
arrive at the following high temperature expression for the internal energy: 



E 

Jr 



6(6 



r(i-S) 
r(i-f) 



t»( T ) ( Y ) 



(6.25) 



Curiously, the T 26 subleading term in (I6.25P vanishes in the unflavored case 6 = 1. Finally, 
from (16.251) we can readily obtain the behavior of the specific heat c v for large T, 



rp 

Af 



(6-l)(3-26) 



r(i-t) 
r(i - !) 



tan ( — ) — 
v 3 ; V T 



(6.26) 



6.3 Speed of sound 

The speed of sound of a thermodynamic system can be obtained from the other thermal 
quantities by the relation 



dP dF (dE\~ l _ s 



(6.27) 



For a conformal system in 2+1 dimensions, as our flavored background, the formula (I6.27P 
yields v 2 s = 1/2. In this section we analyze the effect of the massive flavors introduced by 
the probe in the deviation from this conformal value. With this purpose we will apply (I6.27P 
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to the background plus probe system, i.e., we will substitute in (I6.27P s and c v by Sback + s 
and c Vt back + c v , respectively, where s and c v denote the entropy density and specific heat of 
the probe and (calculated in (16.111) and (I6.17P ) and Sb ac k has been written in (I2.35p . Hence, 
we get 

C-v,back &v 

The specific heat of the background is related to its entropy as c v ^ ac k = 2sb ac k- Moreover, 
in the probe approximation the D6-branes produce a small deviation from the conformal 
behavior. By expanding at first order, we arrive at the following result 



3-26 d 



2 Sb ac k dT 



Af 



cm 



(6.29) 



Taking into account that Af oc T 3 and that m oc T b , we can cast (16. 29j) as 



5vt 



1 3-26 Af 



V '~2 
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d (log c) 
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(6.30) 



Moreover, from (I2.35P and (14.311) . one can verify that he ratio Af /T Sback can be put as 



Af 



A C 



T Sback 4iV&e 



1 A g 2 
4 iV b A ° 



(6.3i; 



where, in the last step, we have used (14. 7p to write the result in terms of the screening function 
cr defined in (I2.30p . Plugging (I6.3ip into (I6.30p . we arrive at the following expression for the 



deviation 5v^, 
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2b) a 2 



16 b 4 



d (log c) 



cm 



(6.32) 



d(\ogm) . 

We plot in Fig.[H]the result of the numerical evaluation of 8v 2 s as a function of the temperature 
for different values of the flavor deformation parameter e. We see that in all cases Sv 2 is 
negative, which implies that the massive flavors reduce the speed of sound. This effect is 
larger as we approach the temperature where the phase transition takes place. Generically, 
Sv 2 decreases as the number of massless flavors (and thus of the deformation parameter e) 
is increased. This is simply a consequence of the fact that we are considering only one D6- 
brane probe and therefore its effect is more and more diluted as Nf — > oo. In order to have 
a better understanding of this behavior let us estimate 5v 2 in the low and high temperature 
regimes. At low T we can use (15. 4p to compute the right-hand side of (16.321) . We get 



Svi 



9 A 



qa 



4 N (2b + 3) b 3 



T 
M 



(T -»■ 0) 



(6.33) 



Thus, we find that the temperature dependence of the deviation from conformality at low T 
(i.e., 5v 2 s ~ T 3 ) does not depend on the number of massless flavors. However, the coefficient 



38 



6v s 2 



1 

0.0 
-0.05 

-0.10 


i i 1 i 


0.4 

"s, 

/ ' 1 
/ ' ' ' 
/ ' ' ' 


i i i 
0.6_ - - - 


i I 

-j 1 j^0j8 


-0.15 




/ 1 ' ' 






-0.20 


\ i 


| i \ 






-0.25 











A — 



Figure 8: We plot the speed of sound versus temperature for black hole (black curves) 
and Minkowski (blue curves) embeddings. The solid (dashed) curves are for the unflavored 
(flavored with e = 1) background. The dotted vertical lines correspond to the locations of 
the first order phase transition. 



multiplying T 3 in ( I6.33P does depend on e and approaches zero as e becomes large. To 
illustrate this fact let us evaluate the leading term on the right-hand-side of (I6.33P when 
e — t- oo for fixed 't Hooft coupling. We get 

5v 2 s ~-l-^(— )\ (T->0, e->oo). (6.34) 



Similarly, for large T we can use (15.1 Op to evaluate (16.321) . 



r 2 3 A q (3 - 2b) a 2 

Sv m — — 

8 N ¥ 



r(| - !) 



2 



tan(^)(-) , (T^oo). (6.35) 



Therefore, 5v 2 vanishes for T — > oo as a power law that depends on the parameter b (5v 2 ~ 
T~ 2b ). In this case the addition of massless flavor produces a faster decrease of 5v 2 with the 
temperature. However, the coefficient of this power law increases with e. Actually, one can 
easily verify from (I6.35P that for large e and T, Sv 2 behaves as 

^ 2 ~-^(99 § , (T^oo, 6^oo). (6.36) 
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7 Summary and conclusions 



In this paper we studied the thermodynamics of flavor D6-branes in the gravity dual of 
Chern-Simons matter theory in three dimensions. The background geometry is a black hole 
of type IIA supergravity with delocalized sources which includes the backreaction due to 
massless flavors. The corresponding metric and forms are just the straightforward T ^ 
generalization of the AdS^ x A4q solution found in [16], in which the deformation due to the 
massless flavors is encoded in the constant squashing factors of the different pieces of the 
metric. We added to this background an additional D6-brane probe, representing a massive 
flavor, and determined its holographically renormalized action, which passed several non- 
trivial tests. We then studied the thermodynamics of this probe in the flavored black hole 
geometry. 

At low temperature the probe brane does not intercept the horizon of the black hole and 
we have a Minkowski embedding while, on the contrary, at high temperatures the brane falls 
into the horizon. At some intermediate temperature the system undergoes a first order phase 
transition which can be interpreted as a meson melting transition. We studied the different 
thermodynamic functions for both types of embeddings, as well as the corresponding phase 
transition. All the results depend both on the temperature and on the flavor deformation 
parameter e oc Nf/k. The dependence on the latter is encoded in the different functions (b, 
q, a, and £) of the background. 

It is important to understand the different scales of our system. We notice that the 
background has only one independent scale, namely the temperature T. A massive flavor 
introduces a new scale in the problem, which is precisely the mass m q of the quarks. This new 
scale is better characterized by the mass gap of the quark-antiquark bound states, defined 
as the mass of the lightest meson at zero temperature. Up to numerical factors this mass 
gap is the quantity M defined in flD.3j) . Notice that M depends on the screening function 
a, which was to be expected since a parameterizes the flavor screening corrections to the 
quark-antiquark Coulomb force. Given these two mass scales of the problem it is very natural 
to consider their ratio. As M /T = m~b , this dimensionless quantity is related to the mass 
parameter m which we used as the independent variable of our thermodynamic functions. 

One of the main targets of the present paper was the study of the impact of the density 
of smeared massless flavors in various observables of the massive probe. The feasibility of 
such investigation relied on the analytic dependence on the deformation parameter e of the 
background. Only a limited number of such observables had thus far been analyzed. One 
such example is the location of the first order phase transition between Minkowski and black 
hole embeddings. As shown in Fig. H] the location increases with T/m q . This has a dual 
interpretation depending on which variable one chooses to keep fixed, either T or m q . The 
increasing behavior is also observed in 3 + 1 dimensions in the particular case of D3-D7-brane 
system which acts as a model for the quark gluon plasma (see the last paper in [20]). Also, 
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the (absolute value of the) condensate |c| increases for any value of m (see Fig. |3]). Notice, 
however, that the relation between c and (O m ) involves the functions b, q, and a (see f)5.2p ). 
The first two functions reach a constant value when e — > oo, while the screening function 
a decreases as 1/ a/I in this limit (see Fig. [I]), which implies that (O m ) — > with infinitely 
many flavors. The c versus m curve enjoys the self-similarity properties in the neighborhood 
of the transition point. We also examined the deviation of the speed of sound 8v 2 s away from 
the conformal result due to the massive probe and found that it decreases as a function of e 

(Fig. ED. 

The work presented here can be continued in several directions. First of all, we could 
study the fluctuations of the probe brane and obtain the meson mass spectrum at non- 
zero temperature. This study would allow to characterize more precisely the meson melting 
transition. Secondly, it is quite natural to analyze the thermodynamics of the D6-brane 
probe at non-zero baryon density and chemical potential which, as in [28] , can be introduced 
by switching on a non-vanishing worldvolume gauge field. A related project would be the 
study of the thermodynamics of the holographic systems introduced in [29], which contain 
self-dual configurations of the worldvolume gauge fields that represent D2-branes dissolved 
in the D6-brane. 

Another possible future direction of the present work could be the addition of a magnetic 
field in order to study the phenomenon of the magnetic catalysis of "chiral symmetry break- 
ing". At weak coupling this has been studied with conventional perturbative field theory 
techniques while at strong coupling a holographic study has been performed using flavored 
A/"=4 Yang-Mills theory [3D]. As a warm-up analysis, the magnetic field could only couple 
with the probe flavor brane while a more elaborate approach would correspond to a coupling 
of the magnetic field with the backreacted flavors of the background (for the similar analysis 
in the D3-D7 case see [3UI32]). 

A combination of the charge density and the magnetic field (with non-vanishing NSNS 
-B-field in the background and/or supplementary internal flux on the worldvolume) would un- 
cover many interesting phenomena with potential applications to condensed matter physics. 
For example, reduced supersymmetry M — 3 — > 1 due to smeared backreacted flavor branes 
may help bypassing the arguments in [J5] and would thereby allow for the study of the 
quantum Hall effect, along the lines of [33]. It would also be important to study how does 
the flavor deformation parameter enter into the physics of the holographic zero sound [31] 
(for T ^ generalization, see |35[I3D]). and, allowing a non- vanishing Chern-Simons term 
as assumed above, the properties of the magneto-roton excitation [37] and the subsequent 
formation of the striped phase away from the quantum Hall phase [35l[38] . 
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A Probe action in isotropic coordinates 



Let us consider the coordinates R and p, defined in (I3.12p . where one should understand that 
u is the isotropic radial coordinate at non-zero temperature introduced in (I4.24p . We will 
first study embeddings of the D6-brane probes that are parameterized by a function R(p). 
The induced metric takes the form 



ds 2 



^ [p 2 + R 2 ] l fl \-Cdt 2 + (dx 1 ) 2 + (dx 2 ) 2 
2a i f 2 



q da 2 + q sin 2 a df3* 



R 2 



b 2 p 2 + R 2 
(dip + cos a d(3 ) 



The determinant of this induced metric is 



V^-det 



7^,3 



L'r 



97 



46 4 



^ q sinap[p 2 + R 2 }-*- 1 ffVl + R'' 



(A.1) 



(A.2) 



where / and / are given in (I4.27P and (14.261) and it is understood that u = yp 2 + R 2 . 
By using these results the Lagrangian density for the DBI part of the probe action can be 
written as 

Cdbi = -Mp[p 2 + R 2 ]^ b ~ l f f VTTR 72 , (A.3) 

where M is the constant which has been defined in (I4.30|) . 

Let us next calculate the WZ term of the action. We first compute the pullback of /C in 
terms of the R and p coordinates for an embedding parameterized by a function R(p). By 
using: 

R — pR' ^ dr 1 / RR' + p 

V ~ b J p 2 + R 2 



de 



p 2 + R 2 



dp 



dp 



(A.4) 
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the pullbacks of e M (ji— 0,1, 2) and e 3 become 



^ /§ [p 2 + i? 2 ]^ 
2§ 



a3 L / RR> + p 
e = — — — dp 



6 f p 2 + R 2 



(A.5) 



while the pullbacks of the other one-forms e a are the same as in (13.261) . Therefore, /C is given 
by 



4 6 4 
A R 



1 - ^ 



fJ p 2 + R 2 



(p& - R) 



d 3 x A dp A S 3 



(A.6) 



where S 3 is the three- form defined in (I2.19p . To evaluate the WZ term we need to compute 
the improving term 8C7. From (14. 14p and (1A.4j) we get 



Li(0) d9 

Moreover, by using the identity, 

1 



r\ Rp 



b (p 2 + R 2 ) 2 



(pR' - R) dp 



(A.7) 



A [p 2 + R 2 f b ff{j~j) =1 



we can insert the unity in f lA.7j) and rewrite this last equation as 



L qrh e-n i (6)de = ^e-*p[p* + R^- 1 ffx 



b 3 



46 4 



x 



L L\ R 

f f)p 2 + R 2 



(pR' - R) dp 



(A.9) 



It is clear by comparing flA.6j) and (IA.9j) that the effect of the improving term L\ is to change 
/// t° ///• By including the zero-point energy term, which is given by 4 6 A J\f, we arrive 
at the following WZ action 



Swz = M I d 



X 



dpp[p 2 + fi 2 ]»-' // 



By adding f1A.3h and (1A. 10h we obtain the total action S in the (p, R) variables 



+ 46A 
(A.10) 



S = -Af d 



dp P [p 2 +R 2 ]^ff [VTTr^- i+(4-ij - j— ^ (pi?' - a) 



/ 



■4 6A 
(A.ll 



43 



The free energy density F written in f)4.32p follows immediately from flA.llj) . Moreover, by 
taking 4 6Ao = 1 the action flA.llj) coincides with the one written in (I4.52p . 

Let us next take the isotropic coordinate u as the independent variable and let us rep- 
resent the configuration of the probe by the function xi u ) = cos8(u). In order to find the 
Lagrangian density in these variables we notice that the induced metric on the worldvolume 
now takes the form 

r2 2 

- Ub f3 

P 



Cdt 2 + (dx 1 ) 2 + (dx 2 ) 2 



L 2 1 - x 2 + u 2 x 2 2 

+ TTT^T ^ dU + 



b 2 u 



+ 



b 2 



qda 2 + q sin 2 a d(3^ 



x 2 )(# 



i-x 2 

cos ad(3y 



where x = dx/du. The determinant of this metric is 



7^,3 



- det g 7 



Ur 



4fe 4 



- q sin a u <> 1 / / ^/l — x 2 + u2 X 2 



Therefore, the DBI term of the Lagrangian density is 

£-dbi 



u» 1 // V 1 ~X 2 + u 2 x 2 



(A.12) 



(A.13) 



(A.14) 



where M is the constant defined in f l4.30p . In order to compute the WZ part we have to 
calculate first the pullback of the calibration form /C. By using 

dr f du 
r 



(16 



X 



du 



X 



f u 



(A.15) 



we find the different e a 's in the (u, x) variables 



„ 3 L f du 



L 



L^q 



sin a 



x 2 dp 



b 



-r^fqda 



sin a>x df3 



L 



X 



du 



— a/1 — x 2 (dip + cos ad(3) 



b Vi^x 2 

From these expressions we can show that the pullback of K, is 



K 



4b 4 J J 



2 / • 

X -UjXX 



d 3 x A dp A S 3 . 



(A.16) 



(A.17) 



Let us next compute the contribution to the action of the term of 8C7 containing the function 
L\. We get 



Li(e)d9 = -fxxdu 



7Z ut ff 
46 



?(L _ L 
/ / 



XX du 



(A.18) 
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and we again see that the effect of adding L\ is equivalent to changing /// by /// in K. 
Taking into account the zero-point energy, and using the value of the constant Ao written 
in (I4.36p . the WZ action becomes 



S 



wz 



Af / d 



x 



duub 1 // 



1 2 / 

1 - X -U^XX 



f 



+ 1 



(A.19) 



Therefore, the total action has the following expression 



S = -Af d 



x 



(A.20) 



B Low temperature (Minkowski embeddings) 



In this appendix we will study, following closely the appendix A. 2 of [23], the Minkowski 
embeddings for high mass (or low temperature), in which the D6-brane probe remains very 
far from the horizon. In this case we have embeddings which are nearly flat (with R(p) 
almost constant). Accordingly, we write R(p) as: 



R{p) = R Q + 8R{p) 



(b.i; 



with Rq being constant and large compared with 5R(p). Let us write the approximate 
Euler- Lagrange equation. First of all, we represent dC/dR' as: 

1 dC 



Af dR> 

where fi(p) and /^(p) are given by 



h(p)d p 6R - f 2 (p) 



(B.2) 



hip) = p[p 2 + R 2 o] 



[p 2 + RI] 



2p 2 R 



(P 2 



Rl) 2 



[p 2 + R 2 o] 



;b.3) 



In flB.2j) and flB.31) we substituted R(p) by Rq after computing the derivative of £ with 
respect to R' . Let us next calculate dC/dR. As in [23], after computing the derivative we 
will neglect the terms with R' and we will substitute R(p) by R . After some rearrangement, 
we get 



J_ dC 

Af dR 



(P 2 + ^) 3 



Rq 



(i + 1) Ri 



[p 2 + Rl] 



(B.4) 



15 



Let us next define a new function h{p) as 

fM = ~ 
This integral can be computed explicitly 
2p 2 Rq l 26 



hip) 



iP 2 + Rl 



2\2 



3 + 26 



(P 2 + ^) 



, 1 dC 
dp MdR- 



(B.5) 



3 + 26 



Rl 



iP 2 + ^) 2+ ^ rI + 



The Euler-Lagrange equation of motion for SR(p) can be integrated once 

h{p)d p 5R + hip) = hip) , 



(B.6) 



(B.7) 



where we have imposed the boundary condition R'(p = 0) = d p 5R(p = 0) = 0. From this 
equation we get 



d p 5R 



h — h 



QRn 



h 3 + 26 P [(p2 + Rtf b - 1] _ 

Then, the asymptotic value of 5R is given by 



1 



1 + ^ 
K 



.2. ii3 



(B.8) 



5R(p — )• oo) 



6i?, 



o 



3 + 26 



F(i? ) 



where F(R ) is, 



^ l - (1 + g 2 ) 1 
6 il + g 2 )~ b 



2b 



Rn 



(B.9) 



(B.10) 



Notice that the integrand in (IB. 101) behaves as g 1 1 for large g and therefore the integral 
only converges if 6 < 3/2 (which is always true because the maximum value of 6 is 5/4). 
At leading order in R we can substitute F(R ) by F(oo) in (IB.9[) . Thus, 



5R(p — > oo 
We find the following value for F(oo): 

F(oo) = -- 



bK ° F(oo) . 



3 + 26 



26 



3-26 



(B.n; 



(B.12) 



where ip(x) = T'(x)/T(x) is the digamma function. For 6 ~ 1 one can represent F(oo) in 
powers of 6 — 1 as 



F(oo) = ~ - log(2) - ™±?* (b - 1) + 



(B.13) 
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The approximate asymptotic value m of R(p) at p = oo can be related to Ro as 



i?o + a (6) i? ( 



) 



where 



a(b) 



6 



Ffoo] 



26 



3 + 26 v y 3 + 26 L 3 — 26 r V6/ r V26 
The relation (IB.14P can be easily inverted at leading order. We find 



R m m — a(b) m l ~ 



(B.14) 
(B.15) 

(B.16) 



In the particular case 6 = 1 the previous formula gives rise to the following relation between 
Ro and m 

Ro~rn - A(3 + 41og2)^ , (6=1). (B.17) 

10 m? 

which should be compared with eq. (A. 12) of [25]. Let us next study the large p dependence 
of R(p). In this limit flB.9|) reduces to 



W R 3T26 



R 



o * P 1 ' 



This equation can be integrated immediately 



5R ~ constant + 



66 



-R 



o ° P 



46 2 -9 

From ( 1B.19P we read the value of the condensate constant c as a function of Ro, 

66 -i-3 
C 46 2 - 9 ' 



(B.18) 



(B.19) 



(B.20) 



In particular, in the unflavored background 6 = 1 we get the following relation between c 
and m: 

c^-?^r, (6 = 1). (B.21) 



5 m 4 



B.l On-shell action 



Let us use the previous results to evaluate the on-shell action for the Lagrangian density 
(lA.lip . To compute the leading order result at high mass (or low temperature) it is enough 
to take R = Rq in the action. Let us express the result in terms of the function Q (m) defined 
in (16.11) . By taking R = R in (16.31) . we get 



G{m) 



2R 2 G 



dpp [p 2 + Rl] 



[ P 2 + RlY b - 1 



(B.22) 
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The integral on the right-hand side of flB.22j) can be easily performed, 

2b 



Q{m) 



1 - 



rJ 



2b + 3'° ' 

At leading order we can take Rq = m in (IB. 231) and we get the estimate (16.181) . 



(B.23) 



C High temperature limit (black hole embeddings) 



Let us now consider the limit of high temperature or, equivalently, low quark mass. Note 
that the D6-brane embedding with x — is an exact solution to the equation of motion. 
In order to study solutions for which \ remains small, we expand the D6-brane action to 
quadratic order in x an d obtain the following equation of motion 



0„ 



_3 

U b 



1 -w) x-uil -u» x 



+ U 1 b ( 1 — U 'b 



l + x-u 1 - uH x 



The general solution of (IC.ip is 



L U~b 1 F(- , 1 



6 3 b e 

-; : u b 

3' 2 3 



„/l b 1 b 6 
CouF — , -; — | — :u b 
\2 ' 3 2 3 



(c.i; 



(C.2) 



where ci and C2 are two constants to be determined. Let us focus on the u = 1 behavior of 
x( u )- in general, we have 



f(oi, /3; a + (3; z\ 



loefl - z) 



r(«)r(/3) 

Therefore, near u = 1 the general solution (10. 2[) behaves as 

1 



as z — > V 



(C.3) 



r(l-l) r , r(j + |) 



r(i-i) 



r(i) 



loefl -ub\ 



(C.4) 



Thus, the solution f lC.2j) is generically singular at the horizon u — 1. To avoid this singularity 
we must impose that the coefficient of the logarithm in (IC.4I) vanishes, which leads to 



a r(i + f)r(i-f) 
r(i)r(l-|) • 

Interestingly, this condition is equivalent to requiring x( u — K 



(C.5) 



0, as in fl£BJ). 
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Let us now look at the behavior at u — oo. To find the asymptotic limit of F(a, (3; 7; z) 
at large z we make use of the following relation 



F(a,f3;r,z) 



r( 7 )r(/j - a) 

r(/3)r( 7 -a) 



-l) a z~ a F(a, a + 1 - 7; a + 1 - (3; -) + 



1 



+ r^r/* m Z ~ P F V> P + 1 " r, P + 1 " «; -) • (C6) 

r<a)r(7-p) 2 

For the particular case of 7 = a + (3 this formula leads to the following asymptotic behavior 
for large z: 



(-1) 



a rpg - a) 
r 2 (/3) 



z" a + (-1)' 



r(a-/3) 

r 2 (c0 



-/3 



(C.7) 



It remains to determine the values of (—1)" and (— I) 13 , which are in general multivalued. 
By comparing with the numerical results when a and /3 are as in (1C.2I) one concludes that 
one should take — 1 = e~ l7T and thus (—1)° = e _;i7ra (and similarly for (— 1)^). Thus, for the 
two hypergeometric functions in (10.2|) . we can write at large u 



~ l F\ - , 1 



uF 



b 3 b 
3' 2 ~ 3 : 



M 6 



-2 3 



■ r(|-|) 1 r(-| + |) ^ 1 



16 1 6 b 
— ; — I — \u b 
2 3 2 3 



r 



1 b 

2 + 3 



r(|-|) c _^i ..r( - 1 + §) 1 



7T 



r 2 (!) 



(C.8) 



Using these equations we see that the coefficient of 1/u in the asymptotic expansion of x{ u ) 
is 



— 1 



r(i-i)r(i-tt , r(i + j)r(i-|) ^ 



(C.9) 



^(l-j) 

The imaginary part of flC.9j) should be zero (otherwise the mass would be complex). This 



condition leads to 



c i = _sin(f) r(§ + f)r-(i-§) 

c 2 



(CIO) 



r(|-|) ' 

which can be shown to be equivalent to flC.5j) by taking z = 6/3 in the reflection formula for 
the Gamma function, namely 



r(z)ra-z) 



7T 



(c.11; 



sm{Trz) 

The only contribution to the real part of flC.9j) comes from the second term, and is given by 



,ixb, r(~ + 1) r(~ - 1) 

cos ( — — — — — — — c 2 . 

3 7T 



(C.12) 
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One can check that the coefficient multiplying c 2 in this last expression is one by using again 
the reflection formula (IC.llj) . Thus, we can identify c 2 with the mass parameter m. Let us 
next study the subleading terms. The coefficient of 

w -3/6+i ig 

ECj-jM-j+j) . £(i+j)I(-i+j) 

— ^ — - - ^-e 3 ci - i— — d/ tVt — c 2 ■ (C.13) 

r 2 (!) 

By requiring the imaginary part of (1C.13P to vanish we get again an expression for Ci/c 2 , 
which can be shown to be equivalent to f lC.5j) by using (IC.llj) . Moreover, the real part of 

(E33} is 

- cos A r d-f) r (-^ + t) Cl , (C.14) 

3 7T 



which can be shown to be equal to c\ by using again ( IC.llj) . Thus, we can identify c\ with 



the condensate c in (I5.7p . From these identifications of C\ and c 2 and their relation (10.5I) it 
follows that, in this low-mass regime, the condensate c is linear in the mass m and is given 
by 

C S3 - U ,I 3 , , ^ 777. , (C.15) 

r(I)r(l-|) 

which is just the expression written in fl5.10p . Let us now find the relation between the value 
of x & t the horizon (xh = x{ u = 1)) an d the mass m. We find that Xh can be simply written 

as 

rd - ±) 

Xh^V^ r( i_l\ rn , (C.16) 

which coincides with (15.91). 



C.l On-shell action 

We can now use the approximate analytic solution found in the previous subsection to 
compute the on-shell action which is needed to evaluate the free energy at high temperature. 
Instead of applying a brute force method let us use the fact that the on-shell quadratic action 
can be computed as the integral of a total derivative (i.e., by taking the appropriate limits 
without the need of performing the integral). Let us consider first the generic case of an 
action of the type 



S 



du 



l-H) 



Fi{u) x 2 + F 2 {u)x 2 + F 3 {u)xx 



(C.17) 



where the F^s are known functions of the radial variable u. The equation of motion derived 
from S is 



d_ 

du 



Fx 



dx~ 
du - 



F 2 X 



1 d 

2 du 



[F 3 x] + t;F 3 



dx 
du 



(C.18) 
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If we rewrite the action as 



du 



d 
du . 



Fix 



dx 
du 



X 



d_ 

. du 



dx\ 
du J 



F 2 X 



+ F 3X 



du 



(C.19) 



then, after using the equation of motion (IC.18I) . the on-shell action can be written as 



G = S 



on— shell 



du 



d 

du 



x ~r 

du 



(C.20) 



Equivalently, we can write Q in terms of boundary values at u = uq and at u = 00 



= *(*5; + 5«* 



(c.21; 



U=UQ 



We will apply this method to compute the function Q(m) defined in (16. ip in the small m 
regime. Thus, we will take Uq — \ and we will identify the function Q of flC.20j) with Q{m). 
By expanding the right-hand side of (16 .4p to quadratic order in x we find, 



FAu) 



F,(u] 



1 3 +1 
- U* + 

2 



1 8+1 1 

2 2 



U b + 



1 1-1 t t 

■K ub ff 



— U b 

2 



- u b 
2 



Notice that FAu 



F 3 (u) 
: F 3 (u 



u» r 



3 

u >• 



_ 3 

u b 



2 . 



(C.22) 



1) = and therefore there is no contribution from the 
horizon to the right-hand side of (1C21|) . Moreover, from the UV asymptotic behavior (15. 7p . 
we get for large u 

Fi^ + If 3X = _^-(3-26)« + ••• , 



du 2 /v 2b 
where the dots represent terms which vanish when u — > 00. Therefore 



lim x ( -^l 



dx 



F3X 



3 -2b 



2b 



cm 



du 2 

and Q(m) can be approximated in this large temperature regime as 

3-26 



Q(m) 



2b 



cm 



(C.23) 



(C.24) 



(C.25) 



which is just the expression used in the main text (eq. (I6.2ip ). Let us rewrite this equation 
in a more explicit way. By using (1C.15P and the reflection formula (1C.11I) . we can write 

2 



3 -2b 



2b 



cm 



3 r(| + |)r(i-|) 
'* r(|)r(i-f) 



m 



r(i-|) 
r(i-i) 



Tibs, 



tan (— ) 



m 



(C.26) 



51 



Therefore, finally we arrive at 



Q(m) 



r(i-|) 



i) 



tan (—J m 
o 



(C.27) 



D Mass and condensate 

In this appendix we study in detail the relation between the parameters m and c, the quark 
mass m q , and the condensate (O m ). The quark mass m q can be obtained by computing the 
Nambu-Goto action of a fundamental string hanging from the boundary to the horizon. The 
relation that is found in this way is 



1 L*r h 



i 



It is easy to write the right-hand side of fID.ip in terms of gauge theory quantities. First 
of all, we recall that = 4-7rT/3. Moreover, the AdS radius L for the flavored background 
is given by L 2 = it a/2A a, where A = N/k is the 't Hooft coupling and a is the screening 
function defined in (I2.30p . By using these equations we can rewrite (ID. 1H as in (15.11) . This 
expression can be inverted, 

3m q l\b f M 



where, in the last step, we introduced the quantity 



M = , 3m « . (D.3) 
2*tx V2X a 



It follows from (ID. 21) that for fixed quark mass m q and 't Hooft coupling A, m depends on 
T as m oc T~ b . 

Let us now turn ourselves to the calculation of the condensate, which can be obtained 
from the derivative of the free energy with respect to the bare quark mass fi q , 

(O m ) = g . (D.4) 

In order to compute the derivative on the right-hand side of (ID.4j) we should find the relation 
between the bare mass \x q and the mass parameter m. Notice that the quark mass written 
in (15.11) contains the screening effects due to quark loops, which should not be included in 
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the bare mass. These effects are encoded in the functions b and a. By taking b = a = 1 in 
(15. ip we switch off the dressing due to dynamical flavors. Accordingly, our prescription for 
the bare mass \i q is 

H q = — - V2A Tm . (D.5) 



By using the chain rule we can relate the derivative with respect to fi q to the derivative with 
respect to m. Actually, it follows from (ID .5^ that 

dm 



dn q 



m 



/'</ 



(D.6) 



m 



Af 



Ik, 



d_ 

dm 



F 
Af 



(D.7) 



and therefore the condensate is given by 

H q om 

where, in the last step, we multiplied and divided by Af, which was defined in (I4.30p and does 
not depend on m. The calculation of the derivative on the right-hand side of (ID.7[) is very 
similar to the one performed at the end of Section [3] in the zero temperature background. 
As in (16.21) . we will represent F/Af by means of the integral Q(m). We will work with the 
(u, x) variables and parameterize these quantities in terms of a density F(u, x, X) as follows 



Jf = G{m) - 1 



du T(u, x, X) ~ 1 



(D.8) 



The explicit expression of fF can be read off from the right-hand side of (16. 4p . Notice that 
Qim) depends on m implicitly through the embedding function xi u )- Indeed, changing m 
is equivalent to modifying the boundary conditions for the embedding, which in turn gives 
rise to a new solution of the equations of motion of the probe. The variation with respect to 
m of F/Afu for a function x that satisfies the equations of motion can be obtained from the 
asymptotic behavior of the derivatives of fF. Following the same steps as in ( I3.45P we arrive 
at 

d / F \ dfF dx 
dm \Af/ dx dm 



(D.9) 



U=UQ 



Let us consider from now on a black hole embedding for which uq = 1 (a similar result can be 
obtained for the Minkowski embeddings by working in the (p,R) variables). The derivative 
of fF appearing in (ID.9P is 



As f(u 



dF ~ 3 +1 
dx 

(see (@2ZD), we have 



X 



\/l - X 2 + u2 X 2 

dF 

dx 

u=l 



2 ,9 

+ uo fx 



(D.10) 



(d.ii; 



53 



and the contribution at the lower limit of (ID. 9}) vanish. In order to evaluate the asymptotic 
value at u — > oo, let us remember that x an d X behave as 



m c 
X + s-r + 

Thus, it follows that 



m , 3 N c 

--2 + ( 1 -J j )— + 

Ub 



u- 



u b+1 x ~ —mu b 1 + (l — -) cu + 
Then, for large u, we get 



3. 
Tr 



u b x ~ mu 6 + cm 4- 



26-3 



CM + 



Taking into account that 
we finally arrive at 



dx , 1 

9m it 

dm \J\fJ dx dm. 



2b -3 



c . 



From (ID.16P we readily get 



dF 

dm 



3 -2b 



■J\f . 



(D.12) 

(D.13) 

(D.14) 
(D.15) 
(D.16) 

(D.17) 



By using flD.17j) to evaluate the right-hand side of flD.7j) . we obtain the relation between 
(O m ) and c that we were looking for 



(O m ) 



3 -2b mj\f 

b fig 



c . 



(D.18) 



Therefore, c is proportional to the condensate (O rn ) as expected. Let us now write this result 
in terms of gauge theory quantities. By using (I4.3ip and (ID. 5}) we can write 



mM 2a n ( „ T ~ 
NT. 



9 b 



(D.19) 



Plugging this result into (ID.18[) and using f)4.7p to eliminate a, we arrive at the formula 
written in f)5.2p for the condensate (O m ). 

Let us now determine the high and low temperature behavior of (O m ). We start by 
considering the behavior for large T (or small m) and fixed m q . In this case c ~ m ~ T~ h 
(see flCH5j) and (jLTTjl ) and thus, 



(O m ) ~ T 



2-6 



(T — >• oo) , 



(D.20) 
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where in the last step we wrote the result in terms of the mass anomalous dimension 7 m = 
b — 1. Thus, the dependence on T of the condensate in this high T regime varies with the 
number of flavors. Actually, it is determined by the mass anomalous dimension 7 m . Clearly 
(O m ) grows linearly with T for the unflavored background, whereas, since ■y m = 1/4 for 
Nf — > oo, the condensate only grows as T 3 / 4 when the number of flavors is very large. 
At low T we found in (15.41) that c behaves as 

c ~ m" 1 -! ~ T b+3 , (T -> 0) . (D.21) 

By using this result in ( 15. 2 j) we conclude that the dependence of (O m ) on the temperature 
for low T is given by 

(O ro > ~ T 5+b , (T — >• 0) . (D.22) 



E Thermal screening 



In this appendix we will analyze two quantities that characterize the screening of quarks in 
the thermal medium of our flavored black hole. We will start by studying the quark- ant iquark 
potential, following the approach of [321110] (see also [H1I12]), m which one considers a 
fundamental string hanging from the UV and penetrating into the bulk. If tq denotes the 
minimal value of the radial coordinate reached by the string, one can show that the quark- 
antiquark distance d on the boundary is given by 



d 



dy 



i vV - y(l - h )^y* - yjl - h Q ) - h c 



where ho denotes 



ho = h(r = r ) 



(E.1) 



(E.2) 



Moreover, we can also compute the energy of the quark-antiquark pair by evaluating the 
on-shell action of the string. This quantity must be regulated by subtracting the energy of 
two straight strings stretching from the UV to the horizon r = r^. The final result of this 
calculation yields 



E, 



<i'i 



2Xa <r Q 



dy 



vV - y(i - ho) - h 



r + r h 



(E.3) 



Notice that, as in the T = case, the screening effect due to the dynamical quarks is given 
by the function a multiplying the square root of the 't Hooft coupling A in (IE.3j) . In order 
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to investigate the departure from the Coulomb behavior due to the finite temperature, let us 
expand flE.lj) and flE.31) in powers of T (or, equivalently of r^) and keep the first non-trivial 
contribution. For the qq distance d we get 



d = 

where J is the following integral: 

dy 



2 V2n 3 / 2 



i\ 2 



r ° r (I) 



J" + 



(E.4) 



J 



y 2 Vv 4 



l + y + y 2 + 2y 3 
y 3 (l + y + y 2 + y 3 



0.093 



(E.5) 



The relation 0E.4|) can be easily inverted to obtain r as a function of rf, 



2 y^TT 3 / 2 1 

r(l) 2 5 



(E.6) 



Similarly, i?™ for low T can be expanded as 



E 



Vx 



a 



r(l) 2 



''0 



J 



where J is defined as the following integral: 

1 



J 



dy 



+ y + y 



V2r h + ■■ 



(E.7) 



i y{^- + y + y 2 + y 3 )\/y~ 4 



0.485 



(E.8) 



Plugging the value of ro given by (1E.6|) in f ]E.7j) we obtain the quark- ant iquark energy as a 
function of d at low temperature, 



E qq (d,T) - E (T) 



2\a 



4vr 3 



7T 



r(!) 



^ + 54 r *- 



T 3 rf 2 



(E.9) 



where £? (^) = t V 2 A a T is the zero-point thermal energy introduced by our regularization 
and we have used J — J = 7r/8. The behavior displayed in (1E.9j) for low T corresponds 
to the one expected for a quark- ant iquark pair screened by a thermal bath. Indeed, we 
readily conclude that the first temperature correction in (1E.9[) makes the force between the 
q and the q less attractive. Actually, as in [HJH2], one can evaluate numerically the exact 
expressions flE.ip and f)E.3[) . One finds that the qq distance d reaches a maximum and, 
actually, the Coulomb-like behavior E q q ~ l/d, valid at low temperatures, ceases to exist at 
high temperatures and the quarks become free due to the thermal screening. 

The second observable measuring the thermal screening of quarks that we will analyze is 
the constituent quark mass M c below and near the critical temperature. According to the 
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standard holographic dictionary, M c is obtained by evaluating the action of a fundamental 
string hanging from the flavor D6-brane down to the horizon. Thus, following a similar 
calculation in [22], let us consider a Minkowski embedding in the (R, p) variables. The 
induced metric on the worldsheet of a fundamental string extended in t, R at p = is 



ds\ 



Rb f f-a dt 2 + 



2l 



dR 2 



Taking into account that in the above metric, / = 1 — R 
following value for the determinant of the induced metric: 



and / = 1 + R- 



2lb 



„ 3 ■ 

FT* 



(E.10) 
we get the 

(E.ll) 



The constituent quark mass M c is minus the action per unit time of the Nambu-Goto action, 



1 

2^ 



det g 2 



1 L 2 r h 



2vr 



2§ 



rt 
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(E.12) 



where we have taken a' = 1 and Ro is the minimum value of the coordinate R reached by 
the brane. When R — > oo, Rq ~ m, and the constituent quark mass M c becomes equal to 
the quark mass m q and we recover the relation (15.11) between m q and m. It follows that the 
ratio between M c and m q is given by 



Mc 



i 

nib 



A? v-n 



(E.13) 



Notice that in flET3]) M c ->■ as we approach the critical solution with Rq = 1. Actually, 
one can show that M c decreases monotonically with T as we approach the temperature of 
the phase transition. This is the expected physical behavior for a free quark in a plasma. 
Moreover, at low temperature we can use the analytic result (IB. 16}) to obtain the first 
screening corrections. We find 



— = 1- ^4-^ + - -=- 



T 



AI 



T \3 



3 V M J 



b 9 J KM J 



(E.14) 



where M has been defined in flD.3|) and a (b) in (IB. 15[) . 



F Critical embeddings 

In this appendix we study the critical behavior of the D6-brane probe in the flavored ABJM 
black hole. Following closely [2l ^l4"3lf4"5] we analyze the brane embeddings near the horizon. 
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It is quite convenient to choose a new system of coordinates in which the induced metric 
near the horizon has the form of a Rindler space. In order to find these coordinates, let us 
expand the radial coordinate r in the near-horizon region as follows 

r = r h + Cz a , (F.l) 

where z is a new coordinate, which is assumed to be small in the near-horizon region, and C 
and a are constants that will be determined by looking at the dr 2 part of the metric ( 12 . 1 \ . 
We notice that the expansion of the blackening factor h(r) is 

h( r ) « — z a + ■■■ . (F.2) 

r h 

Then, the dr 2 part of the metric is 

dr 2 a 2 C „_ 2 , 2 



r 2 h(r) 3r/, 



z a ~ 2 dz 2 + ■■■ . (F.3) 



We fix the constants C and a by requiring that z is a Rindler coordinate and that the 
right-hand side of this equation is just dz 2 , which is achieved when C and a are given by 

a = 2, C = ^ = ttT . (F.4) 

We will also want to explore the region in which 9 is small. Accordingly, let us represent 9 
in terms of a new coordinate y given by, 

V=l, (F.5) 

and approximate sin# 9 = by. Written in the coordinates z and y, the metric takes the 
form 

-(2ttT) 2 z 2 dt 2 + r\ \{dx 1 ) 2 + (dx 2 ) 2 } + dz 2 + dy 2 + 



L 2 



b 2 



da 2 + sin 2 a df3 z 



+ y 2 (dip + cos a df3) 2 + ■•• , (F.6) 



where we keep the leading terms and the dots represent terms that do not contribute to the 
embedding of the D6-brane. Notice from (IF.6I) that only one of the three internal directions 
wrapped by the probe collapses at the tip of the brane y — 0. 

Let us take an embedding characterized by a function y = y(z), which is a appropriate 
to describe a black hole embedding. To compute the induced metric we just substitute 
dz 2 + dy 2 = (1 + y 2 ) dz 2 in the metric (1F.6|) . where the dot denotes derivative with respect 
to z. The DBI Lagrangian density becomes 



C DB i oc zy \Jl + y 2 , (F.7) 
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where the proportionality constant is ItsT r\M r . Let us now consider the WZ term of the 
action. We will expand the WZ Lagrangian by taking into account that z and y are of first 
order and y is of order zero. With these assignments Cdbi is of second order. By inspecting 
(I4.10p we notice that there are two terms to look at. First of all we consider 

3 b 2 r h 

r sin 9 cos 9 d6 « b 2 ' r^y-ydz + — - — z 2 yy dz + ■ ■ • . (F-8) 

The first term in this equation gives a contribution at first order to the Lagrangian which, 
however, is a total derivative and does not contribute to the equation of motion. Thus we 
neglect this contribution. The second term in (1F.8P is of order three and will be neglected. 
The second term in (14.101) to be considered is 

2>b 2 r 3 

r 2 sin 2 6 dr = — — ^ zy 2 dz + ••• , (F.9) 

which is also of order three. Therefore, the total Lagrangian density at leading order is just 
the DBI one written in (]F.7|) . The corresponding equation of motion is 



zyy + {yy - z)(l + y 2 ) = . (F.10) 

The parametrization y = y(z) is appropriate to study the black hole embeddings. In this 
case the differential equation (IF. 10)) must be solved by imposing the following boundary 
conditions: 

y(z = 0)=y h , y(* = 0)=0, (F.ll) 

where y^ characterizes the angle at which the brane reaches the horizon. 

In the case of Minkowski embeddings the appropriate parametrization is z = z(y). The 
brane does not reach the horizon and ends at y = at a point whose distance from the 
horizon is determined by the value of z(y) and y = 0. Therefore, the differential equation to 
integrate is just obtained from (IF. 101) by exchanging y -B- z, namely 

yzz + {zz - y)(l + z 2 ) = , (F.12) 

with the boundary conditions 

z(y = 0) = z h , z(y = 0) = 0. (F.13) 

Notice that the equation of motion (IF.10I) and the Lagrangian (IF.7I) are the same as the 
general expressions in [25] with n = 1, which is consistent with the fact that only one the 
internal directions of the MP 3 cycle wrapped by the brane collapses at the tip. 
Clearly, the system is symmetric under the interchange of y and z 

y z . (F.14) 
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This symmetry exchanges black hole and Minkowski embeddings. This means that for any 
black hole solution y = f(z) there exists a Minkowski solution z = f(y) with the same 
function /. 

The critical solution is the following particular solution of ( IF.lOj) : 

y = z. (F.15) 

In this solution the brane just ends at the horizon z = at the point y = (i.e., with 9 = 0). 
Therefore, the critical solution (1F.15P is the limiting case of both Minkowski and black hole 
embeddings. Notice that it is invariant under the exchange (1F.14|) . 

Let us write the critical solution in terms of the isotropic coordinate u. First, we recall 
the relation between u and r, 



u 



3 

_3_ f r 

2b - 



-) 2 (1 + Vh). (F.16) 



Taking into account that y/h ~ 3z/4, we have 

3 3 

U 2b^l + -z -> u^l + bz. (F.17) 

It follows that R = u cos 9 and p = u sin 9 can be expressed in terms of z and y as 

R^l + bz, p^by. (F.18) 

Then, the critical embedding y = z in the near-horizon region is given by the following linear 
relation between R and p: 

R = 1 + p . (F.19) 

Thus, dR/dp = 1 and the incidence angle of the critical embedding in the (R, p) plane is 
7r/4 for all values of b. 

Let us next analyze the near critical black hole solutions (the corresponding analysis for 
the Minkowski embeddings can be obtained by exchanging y with z in what follows). We 
represent y(z) as, 

V = z + £{z) , (F.20) 

with £(z) being a small function of z. At first order in £, the equation of motion (IF. 101) 

reads, 

z 2 l + 2zi + 2^ = 0. (F.21) 

This equation can be solved by a power law £ = z v where the exponent v satisfies the 
quadratic equation u 2 + v + 2 = 0, whose two solutions are 

v ± = ~±ai, (F.22) 
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with a being 

a = ?! . (F.23) 
The two independent solutions for £ can be taken to be 

12 ±2 

— y=r sin [ a log(T z) J , — -=r cos [alog(T z) J , (F.24) 

where we have introduced the temperature T in order to deal with dimensionless quantities. 
Therefore, we can write, 



y = z + 



T-i 



A sin [alog(T^)] + B cos [alog(Tz)] , (F.25) 



with A and I? being two coefficients. Notice that A = B = for the critical embeddings 
and therefore the coefficients A and B measure the deviation from the solution ( 1F.15|) . Let 
us denote by m* and c* the values of the mass and condensate parameters which correspond 
to the critical embedding, respectively. Clearly, A and B depend on the differences m — m* 
and c — c*. Actually, it was suggested in pE3lf35] that A and B depend linearly on m — m* 
and c — c*. 

The differential equation (1F.10I) satisfies the following property. If y(z) = f(z) is a 
solution of the differential equation (IF. 101) . then y(z) defined as: 

Viz) = ^ , (F.26) 
H 

with fi G R being an arbitrary real number, is also a solution of (IF.lOp with the initial 
condition, 

y h = y( z = 0) = yiL. (F.27) 

n 

Clearly, any two solutions of (1F.10|) with the conditions ( IF. 1 If) are related by this symmetry. 
Thus, we can reconstruct all black hole solutions from a given (fiducial) one. Let us see how 
(1F.26P is realized in the coefficients A and B. First, we define the rotation matrix -M(/i) as 

(cos [alog(/x)] sin [a log(/i)] \ 
(F.28) 
— sin [a log(/i)] cos [a log (/i)] J 

Then, if we denote by A and B the coefficients for the transformed solution y(z), one finds 
that they are related to the initial coefficients A and B by the following combined scaling 

. _ . 

= ^M(fi) f 
^ \ B 
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It is very illustrative to rewrite this result in terms of the z = values of y(z) (yh and yh) 
From flF\27l) it follows that 

(F.30) 



Vh 

H = — . 

Vh 



Moreover, one can check that the matrix M. satisfies, 

M(fi) = M(y h )M-\y h ) . (F.31) 
By using this result in the transformation law flF.29j) we can rewrite this last equation as 

.4 



M(y h 

3 

Vh 



M(y h ) 



B 



Vh 



A 
B 



(F.32) 



where v is a constant vector (it is the same for all the embeddings) . Therefore, the quantity 
on the left-hand side of flF.32j) is the same for all black hole solutions. Let us rewrite this 
property as 

' A 

M-\y h )v . (F.33) 



Vh 



B 



Next, we notice that M.{yh) (and its inverse) is a periodic function of logyh- Actually, it 
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Figure 9: Values of m and c around the critical point for Minkowski embeddings. The solid 
(dashed) curves correspond to e = (e = 10). 

follows from flF.28[) that Ai(yh) does not change when ^- log^ is shifted by one. Therefore, 

_ 3 _ 3 

we get from (1F.33I) that y h 2 A and y h 2 B are periodic functions of f- log yh with period one. 
Since the coefficients A and B are linearly related to m — m* and c — c*, 



logy h ) , —r^ = FA?- logy h ) , (F.34) 

2 V47T / 2 V47T 
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where F m (x) and F c (x) are periodic functions of x with period one. 

A similar result can be found for Minkowski embeddings by exchanging y <H- zq. Then, 
we can write, 



m — m* 

3 
_ 12 



C 

3 
-, - 



Gc (5 iogzh]> ' 



[F.35) 



with G m (x) and G c (x) being periodic in x with unit period. 

Let us recast the previous results in terms of our physical variables. We consider first 
the case of the Minkowski embeddings, which are characterized by the value Rq of the R 
coordinate at p = 0. From the relation ( 1F.18P between R and z, it follows that z and Rq 
are related as 

Ro-l 



By using this result in (IF.35j) . we can write 



m — m* 



= gJ— iogCRo-r 



c c* 



= G c f— log(i2o-l' 

(Ro-iy* 



where the new functions Gm Ax) are defined as 



~ 3 Ol 

G m , c ( x ) = G mfi (x - — log b) . 

Ztx 



(F.36) 



(F.37) 



(F.38) 



It follows from this definition that G m , c (x) are also periodic functions of x with unit period. 
The numerical results for the functions written on the left-hand side of (1F.37|) are plotted in 
Fig. O They confirm this periodicity behavior. 




log(l-^,,) 




log(l-;r/,) 



Figure 10: Values of m and c around the critical point for black hole embeddings. The solid 
(dashed) curves correspond to e = (e = 10). 

Similarly, we can deal with the case of black hole embeddings. In this case the solutions 
are characterized by the value Xh of x = cos 6* at the horizon. For near-critical solutions 



63 



Xh ~ 1 — @h/2 and, by using (1F.5|) . we get that the relation between Xh and is 



J/h = — (1 - Xh) 2 ■ 



Plugging this result in (1F.34|) we find 



m — m* 



:i-xft)* 



V7 



V4tt 
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where the new functions F m ^ c (x) are defined as: 
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(F.39) 



-Xft) ) , (F.40) 



(F.4i; 



Clearly, F m}C (x) are periodic functions of x with period two. The numerical values of (m — 

3 3 _ 

m*)(l — Xh)~ J and (c — c*)(l — x^)~ 3 are displayed in Fig. [TUland agree with the predicted 
periodic behavior. 
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